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ABSTRACT 


The  purpose  of  the  Air  Force  freight  distribution  system  is  to  meet 
spare  parts  demand  requirements  at  minimum  cost.  Budgetary  constraints 
tiaV  suggestfed^hat  total  expected  backorder  level  for  items  at  user  in- 
stallations be  minimized, subject  to  a given  dollar  expenditure  level  for 
inventory  investment.  LOGAIR  (a  dedicated  Air  Force  air  transport  ser- 
vice) is  a major  transportation  sub- system  to  support  spare  parts  delivery 
requirements  of  users  of  high  priority  items.  A two-echelon  inventory 
system  for  spare  parts  delivery  exists  with  centralized,  specialized  in- 

C ^ 

)U 

ventories  at  (the  Air  Logistics  Centers  (ALCs)  and  decentralized,  broad 
profile  inventories  located  at  user  installations.  LOGAIR  provides  trans- 
port with  low  order  and  ship  time  to  reduce  resupply  time  in  the  mainten- 
ance of  inventory  safety  levels  at  user  bases. > A systerns  approach  is 

“ ~Z 2 i.  rn  5 V y 

used  to  formulate  a cost/benefit  modeLwhich  recognizes  the  impact  of 

S&  r 

the  Air  Force  resupply  system  upon'^h^’  total  spare  parts  distribution 
'Sy&teai'in  terms  of  total  inventory  investment  level,  total  system  cost, 
and  backorder  level.  Given  a total  expenditure  level  available  for  alloca- 
tion between  inventory  investment  and  transportation,  the  problem  is  to 
determine  the  optimal  fractional  allocation  to  be  made  to  transportation 
(the  remaining  fraction  to  be  allocated  to  inventory  investment)  such  that 
total  expected  backorder  level  is  minimized. ^A  conceptual  framework 
for  trade-off  analysis  for  minimizing  total  system  cost  in  terms  of  inven- 
tory investment  level  and  resupply  time  level  for  a given  backorder  level 
is  presented.  \ This  conceptual  framework  also  allows  for  the  minimization 
of  backorder  level  in  terms  of  inventory  investment  level  and  resupply 
time  level  for  a given  total  dollar  expenditure  level.  Queueing  models 
exhibiting  various  demand  and  resupply  processes  are  explored  and  com- 
pared to  determine  the  impact  of  inventory  investment  level  and  resupply 
time  level  upon  backorder  level.  Specific  solution  procedures  are  developed 
for  and  are  applied  to  the  trade-off  analyses  mentioned  above. 
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SECTION  1 
INTRODUCTION 


The  purpose  of  the  Air  Force  freight  distribution  system  is  to  meet 
spare  parts  demand  requirements  of  user  installations  at  minimum  cost. 
Budgetary  constraints  have  suggested  that  total  expected  backorder  level 
for  items  at  user  installations  be  minimized  subject  to  a given  dollar  ex- 
penditure level  for  inventory  investment.  LOGAIR  (a  dedicated  Air  Force 
air  transport  service)  is  a major  transportation  sub- system  to  support 
spare  parts  delivery  requirements  to  users  of  high  priority  items.  A 
two-echelon  inventory  support  system  for  spare  parts  delivery  exists  with 
centralized,  specialized  inventories  at  the  Air  Logistics  Centers  (ALCs) 
and  decentralized,  broad  profile  inventories  located  at  user  installations. 
LOGAIR  provides  transport  with  low  order  and  ship  time  to  reduce  resup- 
ply time  in  the  maintenance  of  inventory  safety  levels  at  the  user  bases. 

The  total  systems  problem  is  to  determine  an  optimal  logistics 
support  system  considering  the  cost  trades-off  which  are  possible  among 
inventory,  repair  and  maintenance,  and  transportation  sub- systems  as 
well  as  management  policies  regarding  procurement  and  level  of  service. 
Within  this  framework,  the  transportation  system  is  to  be  developed  to 
utilize  available  transportation  modes  in  servicing  requirements  in  a 
timely,  cost-effective  manner.  Any  mode  which  is  dedicated  to  Air  Force 
needs  (LOGAIR,  truck)  must  be  defined  in  detail  to  include  route  specifi- 
cation, vehicle  schedules,  and  managerial  procedures.  Thus,  a systems 
approach  is  proposed  to  explore  formulation  of  a cost/benefit  model  which 
recognizes  the  impact  of  the  Air  Force  resupply  system  upon  the  total 
spare  parts  distribution  system  in  terms  of  total  inventory  investment 
level,  total  system  cost,  and  backorder  level.  This  approach  recognizes 
that  the  purpose  of  expenditures  for  transportation  is  to  reduce  resupply 
time  in  the  support  of  inventory  policy  designed  to  meet  user  needs. 


Implicit  in  the  approach  taken  by  any  inventory  model  is  the  assump- 
tion of  the  existence  of  a transportation  system  capable  of  delivery  of  an 
order  and  ship  time  (O&ST)  profile  which  makes  up  part  of  the  resupply 
time  required  for  input  to  the  model.  The  Air  Force  freight  transportation 
system  is  supposed  to  provide  timely  delivery  of  inventory  items  to  users 
from  the  supply  depots  and  to  ship  items  from  the  user  installations  back 
to  the  depots  for  repair.  A transportation  system  with  low  average  order 
and  ship  time  levels  should  yield  correspondingly  low  levels  of  required 
inventory  investment  at  a cost  which  is  small  compared  to  any  resultant 
savings  in  inventory  investment.  Another  purpose  of  the  transportation 
system  is  to  minimize  resupply  time  in  the  event  of  a stockout. 

The  problem  then  becomes: 

Given  a total  expenditure  level  available  for  allocation 
between  inventory  investment  and  transportation,  what 
is  the  optimal  fractional  allocation  to  be  made  to  trans- 
portation (the  remaining  fraction  to  be  allocated  to 
inventory  investment)  such  that  the  total  expected  back- 
order level  at  the  user  installations  is  minimized? 

In  this  allocation,  it  is  assumed  that  expenditures  are  to  be  expended  in 
an  optimal  fashion.  Inventory  models  which  can  treat  consumable  as  well 
as  recoverable  items  are  to  optimally  allocate  total  inventory  investment 
among  supply  depots  and  user  installations  subject  to  a given  input  O&cST 
profile.  A transportation  network  and  mode  selection  model  can  be  re- 
quired to  optimally  utilize  a given  transportation  expenditure  level  to 
yield  an  optimum  O&ST  profile  that  minimizes  average  O&ST  level  experi- 
enced within  the  network.  Or,  conversely,  a transportation  network  and 
optimal  mode  selection  model  might  be  designed  to  meet  specific  O&cST 
standards  for  all  shipments  in  such  a fashion  as  to  minimize  total  trans- 
portation cost. 

The  above  approach  ignores  maintenance  and  repair  costs  and  other 
costs  associated  with  other  segments  in  the  resupply  cycle.  Expenditure 
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levels  in  these  areas  are  assumed  to  be  optimal  in  relation  to  the  total 
expenditure  level  allocated  for  inventory  investment  and  transportation. 
In  principle,  a methodology  could  be  developed  to  include  optimal  alloca- 
tion of  resources  to  other  segments  of  the  resupply  cycle. 
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SECTION  2 

GENERAL  MATHEMATICAL  STATEMENT  OF  PROBLEM 


Let  [N  ] be  the  expected  number  of  backorders  for  a given  single 
investment  (recoverable)  item.  [N  ] can  be  written 


[Nb]  = [NB(N,p)] 


(2.  1) 


where 


total  number  of  recoverable  items  including  those  installed, 

X/u.  = XT  = ratio  of  item  demand  rate  to  item  resupply  rate, 
1 v’ 

— = item  demand  rate, 


mean  time  between  demands, 

— = resupply  rate, 

T 

mean  resupply  time. 


The  total  cost  associated  with  maintaining  the  inventory/resupply  system 


can  be  written 


C = C(N,  p ) 


(2.2) 


Suppose  it  is  desired  to  maintain  a given  "safe"  backorder  level  so  that 


[Nb]  - N, 


(2.3) 


where  N ° is  a specified  constant,  and  the  objective  is  to  minimize  sys- 
B 

tern  cost.  Then  the  problem  becomes 


minimize 


C = C(N,p) 


subject  to 


[NJN.p)]  = N(°} 


(2.4) 
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To  facilitate  discussion,  assume  N and  p to  be  continuous  variables;  the 
generalization  to  the  case  vhere  N is  discrete  can  be  achieved  through 
usage  of  first  difference  operations  for  discrete  variables.  Thus, 
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W 


Thus, 


and 


(2.  10) 


(2. 11) 


which  implies 


fCN.p) 


(2. 12) 


If  equation  (2.  12)  and  equation  (2.4)  are  solved  simultaneously,  then  solu- 
tion values  of  (N,  p)  are  obtained  which  minimize  system  cost  C(N,  p).  If, 
on  the  other  hand,  backorder  level  is  to  be  minimized  subject  to  a given 
total  dollar  investment  level  so  that 

C(N,p)  = CQ  (2.13) 


then  equation  (2.  12)  and  equation  (2.  13)  are  solved  simultaneously  to  ob- 
tain the  desired  values  of  (N,  P ). 


In  order  to  determine  the  optimal  (minimum)  system  expenditure 

level  to  be  allocated  between  inventory  investment  level  (N)  and  resupply 

time  level  (T  = pT,  ),  the  following  functions  must  be  determined: 

(0.  A. 

[nb]  - [hb(N,P>]  > 0 ’ (2-  14) 


gl(N,p) 


(2.  15) 


g2CN,P) 


> 0 


(2. 16) 
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MB’S)  ■ S)w 

C - C(N,p) 

8»W.P)  > 

8iW,0)  ■(£)„  < 

8s(N,p)  ■ (Hje  > 

§►'-  The  determination  of  relations  (2.  14)  through  (2.  17)  (STUDY  NO.  1)  is  a 

K»:7 

study  of  the  impact  of  inventory  investment  level  (N)  and  resupply  time 
level  (T  = pT^)  upon  backorder  level.  The  determination  of  relations 
(2.  18)  through  (2.  21)  (STUDY  NO.  2)  is  a study  of  the  impact  of  inventory 
investment  level  and  resupply  time  level  upon  total  system  cost.  The 
results  of  STUDY  NO.  1 and  STUDY  NO.  2 can  be  used  to  determine 
optimum  allocation  of  resources  between  inventory  investment  level  and 
resupply  time  level  as  outlined  above. 

Note  that  the  resupply  time  T is  given  by 

r = (repair  time)  + (transportation  time)  + (administrative  time  for 
^ processing  shipments). 

If  the  subject  of  interest  is  the  impact  of  transportation  time  upon  total 
system  cost  and  required  inventory  investment  level,  then  the  cost  function 


> 0 


(2.  17) 


> 0 , 


0 , 


0 , 


0 


(2.  18) 


(2.  19) 


(2.20) 


(2.  21) 
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C(N,  p)  need  only  include  costs  associated  with  inventory  investment  level 

and  transportation  ship  time  level,  where  a differential  change  in 

p (6p  = \6t  ) is  considered  to  be  due  to  a differential  change  in  transporta- 

M. 

tion  ship  time  level. 

The  above  problem  formulation  is  generalized  to  include  base  level 
repair  and  depot  level  repair  (which  results  in  the  inclusion  of  transporta- 
tion time  as  part  of  the  resupply  time),  the  interaction  of  more  than  one 
base  with  a single  depot,  and  can  be  modified  to  include  the  interactions 
due  to  movement  of  many  investment  and  consumable  items  in  the  trans- 
portation pipelines  among  user  bases  and  supply  depots.  If  desired,  the 
impact  of  administrative  processing  time  and/or  of  repair  time  upon  total 
system  cost  in  relation  to  resulting  required  inventory  investment  level 
can  also  be  determined  using  the  conceptual  procedure  outlined  above. 
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SECTION  3 

DEMAND  AND  RESUPPLY  PROCESSES  IN  INVENTORY  MODELING 


When  dealing  with  failures  of  units  in  operation  it  is  reasonable  to 
assume  that 

1.  The  probability  for  a given  unit  in  operation  to  fail  in  any  parti- 
cular operational  time  interval  of  arbitrary  size  is  the  same  for 
all  identical  operational  units; 

2.  The  fact  that  a given  unit  in  operation  has  failed  in  a given  opera- 
tional time  interval  does  not  affect  the  probability  that  other 
identical  units  in  operation  may  fail  in  the  same  operational 
time  interval  (independence  of  failures); 

3.  The  probability  for  a given  unit  in  operation  to  fail  during  a 
given  operational  time  interval  is  the  same  for  all  operational 
time  intervals  of  equal  arbitrary  size  (the  mean  operational 
time  to  failure  is  long  compared  to  the  total  operational 
time  period  of  observation)  no  matter  when  in  time  the  failure 
takes  place  (time  independence  of  failure  rate). 

The  above  conditions  are  necessary  and  sufficient  to  define  the  probability 
for  failure  of  a given  unit  in  operation  during  a given  finite  operational 
time  interval  t.  Thus  for  small  At,  let 


XAt  = probability  for  failure  of  an  item  in  operation  during 
the  operational  time  interval  At  , 


(3.  1) 


1 - XAt  = probability  of  a given  unit  in  operation  for  not  failing 
during  the  operational  time  interval  At 


(3.2) 


If  P(t)  is  the  probability  that  a given  unit  in  operation  does  not  fail  during 
a finite  operational  time  interval  t,  then 


P(t  + At)  = (1  - XAt)  P(t) 


(3.  3) 
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so  that 


dP(t) 

dt 


- lim 
At  -*•  0 


P(t  + At)  - P(t) 
At 


- - XP(t)  ; 


(3.4) 


dP(t) 

dt 


(3.5) 


(3.6) 


The  above  results  apply  to  a single  unit  in  operation  and  are  independent 
of  the  number  of  units  in  operation  in  the  system.  If,  in  addition,  one 
requires  that 

4.  The  total  number  of  units  in  operation  and  the  total  number  of 
equal  arbitrarily  chosen  operational  time  intervals  are  large 
(making  statistical  averages  significant), 
then  the  failure  process  is  Poisson.  If  the  probability  for  observing  m 
failures  in  an  operational  time  t is  given  by  P^t),  then 


? (t  + At)  - Pa(t)(l  - XAt)  + P^UOXAt  , « > a > l ; (3.  7) 


Po(t+At)  - P0(t)(l-XAt) 


Thus  when  At  -*  0,  one  obtains 


a ■ 0 


(3.8) 


- »,<«> 


m - 0 ; (3.  9) 


m 


- - XP  (t)  + XPm  , (t) 

a m—  l 


> m > 1 . (3.  10) 
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Application  of  standard  techniques  obtains 

P (t)  - . 0<m<«  . (3.11) 

m mi 

It  is  observed  that  the  Poisson  process  only  applies  if  the  number 
of  units  in  operation  is  very  large,  or  if  as  soon  as  a unit  fails,  it  is 
immediately  replaced  (resupply  time  = 0)  by  an  operational  unit.  If  no 
backorders  occur  (i.e.,  there  always  exists  sufficient  servicable  units 
in  stock  available  for  necessary  immediate  replacement  of  failed  units), 
then  the  Poisson  failure  process  applies;  when  backorders  are  allowed  to 
exist  for  time  intervals  of  significant  duration,  or  if  replacement  time  is 
significant,  then  the  Poisson  failure  process  does  not  apply.  Thus  while 
exponential  operational  failure  times  may  be  assumed,  in  general  a 
Poisson  failure  process  cannot  be  assumed  without  careful  investigation 
into  the  nature  of  actual  resupply  times. 

Assume  that  probability  to  complete  a service  on  a given  unit  satis- 
fies the  following  conditions. 

1.  The  probability  for  a given  unit  in  service  to  complete  service 
in  a particular  time  interval  is  the  same  for  all  identical  units 
in  service  being  served  by  identical  servers; 

2.  The  fact  that  a given  unit  in  service  has  completed  service  in  a 
given  time  interval  does  not  affect  the  probability  that  other 
units  in  service  in  other  servers  may  complete  service  in  the 
same  time  interval  (independence  of  service  completion  times); 

3.  The  probability  for  a given  unit  in  service  to  complete  service 
during  a given  time  interval  is  the  same  for  all  time  intervals 
of  equal  size  (the  mean  time  to  complete  service  is  long  com- 
pared to  the  total  period  of  observation)  no  matter  when  the 
service  completion  takes  place  (time  independence  of  service 
rate); 

4.  The  number  of  units  in  service  does  not  affect  the  probability 
for  completion  of  service  for  an  item  in  service. 
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The  above  conditions  are  necessary  and  sufficient  to  define  the  probability 
for  completion  of  service  for  an  item  in  service  during  an  arbitrary  given 
time  interval  t.  Thus,  for  small  At  let 


uAt  = probability  for  completing  service  for  an  item  in 
service  during  At  , 


(3.  12) 


1 - p, At  = probability  for  not  completing  service  for  an  item  13j 

in  service  during  At. 

An  argument  similar  to  that  given  above  for  failures  obtains  exponentially 
distributed  service  times.  (Note  that  if  the  number  of  units  awaiting  ser- 
vice for  a given  single  server  were  always  non- zero,  then  the  number 
completing  service  in  an  arbitrary  time  interval  would  be  Poisson  distri- 
buted. ) 

In  applying  servers  with  exponential  service  times  in  modeling  the 
resupply  process,  the  nature  of  the  server  configuration  is  extremely 
important.  A few  possible  models  for  resupply  are  given  below. 

If  a single  level  resupply  system  function  is  assumed,  then  resupply 
can  be  described  by  a singly  indexed  resupply  state  so  that 


where 


n p , 
max 


0 s n s n 


N 2 n * n 


1 £ n ^ N 
max 


(3.  14) 


(3.  15) 
(3.  16) 


resupply  rate  when  n units  are  in  resupply, 

constant  = resupply  rate  when  a single  unit  is  in  resupply, 

number  of  identical  parallel  servers  assumed, 

total  number  of  units  in  the  system  (including  installed, 
stock,  and  resupply  items). 


The  single  index  for  resupply  is  n,  the  number  of  units  in  resupply. 

Note  that  when  n =1,  then  a single  server  model  is  assumed,  and 

max 

when  n > 1,  then  a multiple  server  model  is  assumed.  If  n is 
max  max 

chosen  so  that  n = N,  then  no  service  queue  is  formed;  if  n < N, 
max  max 

then  service  queues  do  form. 


For  a dual  level  resupply  system  (two  level  resupply;  e.  g.  , base 
level  and  depot  level)  two  server  models  may  be  postulated.  The  first 
is  described  by  a resupply  state  with  a single  index  so  that 


where 

f 

1 

f 

2 


H-l 

^2 


pn=  np  , 

0 £ n £ n 

max 

9 

(3.  17) 

^n  nmax*1' 

n £ n £ N 

max 

9 

(3.  18) 

1 * n £ N ; 

max 

(3.  19) 

p = f p + f p ; 

11  2 2 

(3.20) 

= probability  for  level  1 resupply  (e.  g.  , base  level  resupply), 

= probability  for  level  2 resupply  (e.  g.  , depot  level  resupply), 

= 1, 

= resupply  rate  for  one  unit  in  level  1 resupply, 

= resupply  rate  for  one  unit  in  level  2 resupply. 


Again  the  single  index  for  resupply  is  n,  the  total  number  of  units  in  re- 
supply. The  form  for  the  above  dual  level  resupply  is  identical  to  that 
for  single  level  resupply  except  that  the  effective  service  rate  y is  given 
in  terms  of  a weighted  combination  of  yx  and  y2  . 

A second  approach  to  model  a dual  level  resupply  system  is  to  use 
a two  index  resupply  state.  The  description  of  resupply  rates  is  given  by 


0 < a < a 

max 


i 


oaxp 


max 


a < N 


(3.21) 
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9 


Up2)  - 6u, 


o < e < s. 


max 


.(O 


emaxU2  » 


8 < 8 < N 

max 


U = fjUj  + f2u2  , f,  + f2  - 1 j 


(3.21) 


where 


a 


resupply  rate  for  level  1 resupply  when  ci  units  are  in 
resupply  at  level  1, 

resupply  rate  for  level  2 resupply  when  8 units  are  in 
resupply  at  level  2, 

number  of  units  in  resupply  at  level  1, 


8 = number  of  units  in  resupply  at  level  2, 

a + 8 = n = total  number  of  units  in  resupply,  0 £ a + p = n £ N, 


a 

max 


max 


^2 

N 


probability  for  level  1 resupply, 
probability  for  level  2 resupply, 

assigned  number  of  identical  parallel  servers  at  level  1, 

1 * a s N, 
max 

assigned  number  of  identical  parallel  servers  at  level  2, 

Up  s N, 
max 

resupply  rate  at  level  1 when  one  unit  is  in  resupply  at  level  1, 

resupply  rate  at  level  2 when  one  unit  is  in  resupply  at  level  2, 

total  number  of  units  in  the  system  (including  installed  items 
as  well  as  items  in  resupply  and  in  stock). 


The  two  index  resupply  state  is  labeled  by  (a,  9),  the  number  of  resupply 
items  at  level  1 and  the  number  of  resupply  items  at  level  2,  respectively. 


In  the  single  index  resupply  state  models  given  above,  the  arrival 
rate  for  resupply  (failure  rate)  can  be  modeled  by 


Xffl  - mQX  , m0  < m < N 


(3.  22) 
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described  by 


where 


.(1) 

m 

.<2) 

m 


m 


= failure  rate  when  one  unit  is  installed, 

= resupply  arrival  rate  for  level  1 when  m units  are  installed, 

= resupply  arrival  rate  for  level  2 when  m units  are  installed, 

= probability  for  level  1 resupply, 

= probability  for  level  2 resupply, 

= number  of  installed  units. 


= minimum  number  of  installed  units  resulting  in  meeting 
total  operational  requirements. 

(l) 

Define  71  to  be  the  resupply  arrival  rate  at  level  1 when  the  total  num- 

” (2) 

ber  in  resupply  is  n and  Tj  to  be  the  resupply  arrival  rate  at  level  2 
when  the  total  number  in  resupply  is  n,  0 < n = a + 9 < N.  Then 


• 


- fjV 

t 

0 < n < N _-  m„  ; 

(3.  27) 

- f2o0X 

9 

0 < n < N - mQ  ; 

(3.  28) 

■ 

Ml 

*2 

i 

n)X  , 

N — m + 1 < n < N ; 

0 

(3. 29) 

- f2(N  - 

n)X  , 

(3.  30) 

The  state  transition  probabilities  for  a singly  indexed  resupply 

state  are  given  for  n = N and  4 = nn,  0 < n < N (to  first  order  in  At)  by 

max  n 


P#(t+At)  “ P,(t)[l  - m9AAt]  + Pj(t)[yAt] 


n - 0 ; (3.3  1) 
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PBU>[1  - m^t-nuAt]  +PBH(t)[(n  + l)u4t]  + (t) , 


pQ(t  + At)  - 


0 < a < N - m0  ; (3.  32) 

PQ(t  + At)  - Pn(t)[l  - (N  - n)XAt  - nuAt]  + PQ+1(t)[(n  + l)yAt] 

I 

+ Pn_x  Ct)[(N  - n + l)AAt]  , 

N-m5+l<n<N  ; (3.  33) 

PN(t  + At)  - PN(t)[l  - NuAt]  + Pjj_j(t)[Mt]  , n-N  . (3.34) 

The  differential  equations  which  result  when  it  -*  0 and  the  corresponding 
steady- state  solution  are  given  in  the  section  entitled  "Conservative 
Parallel  Server  Inventory  Model.  " 

When  y = y for  n > 1 and  y =0  for  n = 0,  the  state  transition  proba- 
n n r 

bilities  for  a singly  indexed  resupply  state  are  given  (to  first  order  in  it)  by 

P„  (t  + At)  - P0  (t)[l  - m0XAt]  + Pj  (t)[yAt]  , n-0;  (3.35) 

PQ(t  + At)  - PQ(t)[l  - m9AAt  - yAt]  + PQ+1  (t) [yAt]  + Pn_x  (t)  [mflXAt]  , 

0 < n < N - m0  • (3.  36) 

Pn(t  + At)  - PQ(t)[l  - (N  - n)XAt  - yAt]  + Pq+1 (t) [yAt] 

+ PQ_j (t) [ (N  - n + l)XAt]  , 

N - m0  + 1 < a < N ; (3.  37) 
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PN(t  + At)  - PNCt)[l  - yAt]  + PN-1(t)[XAt]  , n-N.  (3.38) 

The  differential  equations  which  result  when  At  - 0 and  the  corresponding 
steady- state  solution  are  given  in  the  section  entitled  "Conservative  Single 
Server  Inventory  Model.  " 

The  two  index  resupply  state  model  is  described  (to  first  order  in 

At ) by 

P0>B(t  + it)  - P0ig(t)  [1  it  - it  - „0>  it  - Uj<*>it] 
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The  two  index  resupply  state  system  described  above  is  one  that 
most  closely  describes  the  two-echelon  inventory  system  currently  adopted 
by  Air  Force  management  [88,89,  104]  and  is  of  major  usefulness  in  the 
trade-off  analysis  described  earlier  in  this  document.  Below  is  diagram- 
med the  conservative  two  level  parallel  server  inventory  model,  and  a 
mathematical  treatment  showing  the  equivalence  of  this  model  to  the 
singly  indexed  parallel  server  model  follows. 


CONSERVATIVE  TWO  LEVEL  PARALLEL  SERVER 
INVENTORY  MODEL 


Define 


.CO  i £ 


a-o 


P (t) 
a,n-a 


(3.41) 


and  observe  that 


Pa,n-c.(t>  - (°K<t) 


°<a<n,  0 < n < N, 


(3.42) 


, ..  m.LJ-P.  i.ill  I 


if- 


where  P ^ = the  probability  that  there  are  n units  in  resupply  at  time  t. 


n = ct  + 8,  andy  - • 


Further, 


l(l)  + r/2)  * m0A  = A , 0 < n < N-m0  , 


n(l)  + n(2)  = (N-n) A = A , N-m0  + 1 < n < N 


(3.43) 


Therefore, 


n(1)  + ni2)  = K 


0 < n < N , 


(3.44) 


where  A = resupply  arrival  rate  when  a total  of  n items  are  in  resupply, 
n 

Also  note  that 


,(0 


aux  , 0 < a < n , 0 < n<N  , 


(3.45) 


y(2)  =»  (n-a)y  , 0 < a < n , 0 < n<N  , 


n-a 


when  a = N = S (infinite  number  of  parallel  servers  in  resupply  at 
max  max 

level  1 and  at  level  2).  Observe  that 


(3.46) 


<"+ou1f1(”)f?frc,i,n+1<t> 
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n 

*•  ~ Ua+iPa+i  ,n-a^^  = (n+OUjf  XPQ+1  (t) 


Similarly, 


and 


p (t) 

n-a+i  a,n-a+iv  1 


(n-a+.)U2£,£^+‘(n+‘j 


<n+Ou2f  J2)  ! 


a»o 


*,n-«.Pa,n-a+i<t)  ’ 


.(>) 


“a  pa,8(t)  ■ °<‘.pc,,a-a<t> 


,a-n-a/n\„  , . 
ayififi  (a  Pn(t) 


:)r*rv> 


£ ua‘>Pa,n-a(t)  " "M.V') 


a=o 


y(2)p  (t) 

n-a  a,n-a 


(n-a)y2r 


e /n-il  ,a,n-a-i„  , . 
ny2f2(  a / f if2  Pn(t 


<7  y(2?P  (t)  - nu  f P (t) 

n-a  a, n-a  2 2 n 


(3.47) 


(3.48) 


(3.49) 


(3.50) 


(3.  51) 


(3.  52) 


(3.  53) 
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Continuing, 


<£X-iPa-i,n-a(t> 


(t) 


(3.54) 


£ n(2)P 

^ n-i  ct,n-ct-i 


(3.55) 


so  that 


^ r/^p  (t)  + n^p  (t)l 

",  n-i  o-i,n-cC  V-i  a,n-a-iigJ 


X P (t) 

n- i n- i 


(3.56) 


Therefore,  summing  equation  (3.  40)  over  a with  S = n-a  obtains 


3P  (t) 
a 


X P 00 

n n 


nuPQ(t)  + (n+OuP^U)  + Xn_jPn-i(t)  ’ 


0 < n < N , 


(3.57) 


where  P ' - 0 for  n < 0 and  for  n > N. 

n 


Thus,  the  two  index  resupply  state  model  described  by  equation  (3.40)  re- 
duces to  the  one  index  resupply  state  model  given  by  equation  (3.  57),  and 
is  represented  by  the  model  described  in  the  section  entitled  "Conservative 
Parallel  Server  Inventory  Model"  and  by  equations  (3.  31)  through  (3.34) 


in  this  section. 
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SECTION  4 

CONSERVATIVE  PARALLEL  SERVER  INVENTORY  MODEL 

Consider  the  system  for  resupply  of  recoverable  items  that  is  dia- 
grammed below. 


The  parameters  are  given  as  follows: 

N = total  number  of  items  in  the  system, 
n = number  of  items  in  resupply, 
m = number  of  items  installed, 
s = number  of  items  in  stock, 

X = failure  rate  for  a single  item  while  installed, 

X = failure  rate  of  installed  items  when  number  of  items  installed 

m 

is  m, 

m = minimum  number  of  serviceable  installed  items  required  to 
9 

accomplish  mission  objectives  (m>m8  implies  that  service 
performance  level  is  being  met;  m<  mQ  implies  that  the  in- 
stalled serviceable  items  are  working  at  full  capacity  and 
are  not  fully  meeting  mission  requirements), 
m = minimum  number  of  installed  serviceable  items  resulting  in 


no  backorder  ( m < implies  that  a backorder  exists,  N > 

m > m > 1)» 

1 l ~ ° “ 

p,  = — = resupply  rate  for  a single  item  in  resupply, 

y 

t = mean  resupply  time  for  one  item  in  resupply, 

M- 

o = Xt  = ratio  of  failure  rate  to  resupply  rate  for  one  item  in 

M. 

resupply, 

p,  = resupply  rate  when  n items  are  in  resupply, 
n 


For  a conservative  system  observe  that 


N = m + s + n , 


(4.  1) 


and  require  that 


X - m X , 
mo* 


0 < n < N - m„  ; 


“ (N-n)X  , N-m0+l  < n < N ; 


Un  - ny  , 


0 < n < N 


(4.2) 


(4.3) 


(4.4) 


Thus,  it  is  assumed  that  there  are  N parallel  servers  in  resupply.  The 
state  transition  probabilities  are  given  by 


3P0(N)  (t) 


-m0XP  (t)  + yPjCt) 


n • 0 ; 


(4.  5) 


3P00  Ct) 

M 


- -a  XP  Ct)  + Cn  + l)yP  +1(t)  - nyPnCt)  + WiCt)’  (4'  6) 


0 < n < N - m ; 

o 


4-2 


(4.  7) 


3p£N) (t) 


-(N-n)APn(t)  + (a  + l)yPQ+l(t)  - nyPn(t)  + (N-irH)  (t) , 


N-m0  + l<n<N; 


^N)Ct) 


- NyP^Ct)  + XPN-i(t)  , n - N 


(4.8) 


In  the  above,  P ' Mt)  is  the  probability  that  n items  are  in  resupply  at 
n 

time  t.  Using  standard  techniques,  one  obtains  the  following  stead- state 
solution  (t  = <*>). 


(m0p) 


0 < n < N - m ; (4.  9) 


(mop) 


— V 0n-(S-mo)  (N) 

n!(N  - n) ! P 


N—  a+l<n<N;  (4*  10) 


• (B)P)  ntT,t+l 


(4. 11) 


Equation  (4.  9)  and  equation  (4.  10)  can  be  written 


cto  <V>  on 

N-I  (N  - x) ! » 


P®  - 0v>"  ‘>  75-A).  X, 


«.  < x < N , 


(4.  12) 


0 < x < m - 1 , (4.  13] 


— . 


and  equation  (4.  11)  can  be  written 


p(N) 


N-m 

£ 


n-Q 


0 (m0p) 
a! 


n 


+ (m0p)N 


x-o 


(4. 14) 


The  resupply  states  corresponding  toO<n<N-m  (states  having  full  capa 

bility  to  fullfill  mission  objectives)  would  be  described  by  a Poisson  distri- 
(N) 

bution  if  PQ  were  given  by 


p(N) 


-m  p m Xt,, 

o o V 

e = e 


(4.  15) 


However,  the  Poisson  process  applies  only  if  N - mo  "*  “*  a condition  which 
implies  that  the  probability  to  fail  to  meet  mission  objectives  is  zero.  Also 
observe  that  the  resupply  states  for  N>n>N-m„  +1  (those  states  having 
less  than  full  capability  to  meet  mission  objectives)  are  not  described  by  a 
Poisson  distribution,  resulting  from  the  fact  that  N - mo  is  finite. 


As  observed  above,  the  number  of  backorders  is  zero  for  N > x > mt 

(0  < n < N - nij  ),  and  the  number  of  backorders  is  n1  - x = mx  -(N  - n)  for 

0 < x < mx  - 1 (N  - mj  + 1 < n < N).  Thus,  the  expected  number  of  back- 
(N) 

orders  [N  ] is  given  by 
B 


[n£N)]  - 2 <®1  - *>  PN- 


00 


(4. 16) 


x » 


C»f>] 


(■.p) 


(»,P)HpP)  ( 2?  (W7ix-  ♦ ^ TSHT 


(4.  17) 
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Inspection  of  equation  (4.  12)  and  equation  (4.  13)  obtains 


and  therefore 


[N«"l  • [n<N»]  P^'1'  - (N-m, ) P<N> 


(4.24) 


Thus, 


3[Nf»] 

9P 


KN-U1 


[Nf] 


(4.25) 


Define 


A[<>] 


and 


AP 


(N) 


[Nf+1»]  - [N™] 


p (N+l)  _ p (N) 

s s 


(4.  26) 


(4.27) 


to  obtain 


ap 


1^  ifN'N-1’]  , 

P B 


(4. 28) 


and 


9P 


(N) 


3p 


("<N)1  i p(N-l ) 
P S 


(4.29) 


Also  observe  that  equation  (4.  12)  and  equation  (4.  13)  yield 


3P(N) 
1 I 3 N-x 


N - x + _1 


3P 


(N) 


N-x 


N 


P V3P  / u 


0 < x < N j 

(4.30) 


3P 


(N)\ 


P 


(4.31) 
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.... — 


It  is  observed  that  when  m0  = 1,  then  = \ for  all  m >0,  corre- 
sponding to  the  situation  where  the  failure  rate  is  constant  and  independent 
of  the  number  of  installed  serviceable  units  (except  for  m = 0).  Thus, 


».°°  - £ *.<N>  • 


m - 1,  0<  a<  N ; 


(4.32) 


-dl  nl 


(4.33) 


']  - £ (mrx)  T^r 


(4.34) 


Note  that 


11m  p' 


-P  -Xtu 

e ■ e H 


(4.35) 


11m  F«  . £.~\ 

1}  00  11  11  J 


lia  [N<N)]  - lla  £1  <®r  x>  of^loT  ?o(N) 

N -*■  “ N -►  ® x-o 


(4.36) 


0 for  m ■ 1 and 


p < 1 . (4.  37) 


Thus,  a Poisson  demand  process  occurs  for  large  N,  and  for  large  N and 
p < 1 the  expected  number  of  backorders  is  zero. 


When  m9  = nij  = 1,  then 

PuWi  . s!  P(N)  . p(N) 
£NB  J Nl  P"  PN 


m0  - - i ; 


(4.38) 
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When  N = m0=m1  = p = l,  then 


3[^°] 


(4.46) 


C»,C2>]  - C*^]  - 


(4.47) 


‘ T • 


(4.48) 


A case  of  particular  interest  occurs  when  N = m . =m1  . This  is  the 
case  when  there  are  no  extra  spares  and  when  one  or  more  units  in  re- 
supply results  in  system  performance  degradation.  When  N = = nij , 


pCN)  _ 1 

° Cl  + p)N 


(1  + p)  n!(N  - n) ! 


m.  - N ; 


m0  - N ; 


(4.49) 


(4. 50) 


P?5]-  pOT]  - 


1 + p 


®1  " m0  " N i 


(4.51) 


When  there  is  one  extra  spare,  then 


p(N+0  . Ntl_. 

# 1 + H(1  + p)1 


m0  - N ; 


(4. 52) 


p(N+l)  . N(N  + D- 1 


1 + N (1  + p) 


N+i  nl(H  + 1 - n) i 


, 1 < n < N + 1,  m0  - N ; (4.53) 


(N+l),  . h{  1 - (N  - q)(1  + P)N) 


IK  ' ] 


1 + N(1  + p) 


N 


®i  - »o  ■ N ; 


3p 


H 


[*<N)]  - [N^]  « N 


r 


(N  (1  + p)‘ 


(1  + p)N+1  - 1 


mi  ■ ®o  ■ N ; 


L (1  + P)[l  + N(1  + p)N+l]J 


, H - m0  - ax. 


Thus, 


6[N<N)] 


N 


(1  + P) 


— for  N - constant  , N - o0  ■ ax  ; 


- 


- N 


(1  * o)1*1  - 1 


L a + p)[i  + nu  + p)N+‘]j 


for  p « constant. 


N - m0  ■ nl 
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SECTION  5 
POISSON  MODEL 


K 


In  the  Poisson  Model  the  demand  process  described  in  the  section 
entitled  "Conservative  Parallel  Server  Inventory  Model"  is  replaced  by  a 
Poisson  process  as  was  done  by  Sherbrooke  [104]  in  the  development  of 
"Metric:  A Multi-Echelon  Technique  for  Recoverable  Item  Control"  and 
by  Muckstadt  [88,89]  in  the  development  of  "Mod-Metric:  A Multi-Item, 
Multi-Echelon,  Multi-Indenture  Inventory  Model.  " The  following  for- 
mulae apply  when  this  approach  is  taken. 


p = e'm°p  (m°P) 

n 


n 


0 < n < 00 


(5.  1) 


n ! 


KN|] 


00 

i 


(n-N+m  ) P 

l n 


(5.  2) 


n = N-m1  +1 
(N) 


[NjP  i 4(N^N)] 


-P 


(N) 


(5.3) 


3[N<N>] 


6P 

- = m P ' ; 

o s 

(5.4) 

3P<N> 

s 

3 P 

moPN-mj+l 

(5.5) 

ap(N) 

s 

, p (N+1>  . p(N)  . 

s s 

_p 

N-mj+1 

; (5.6) 

tNBN>> 

r (N).  (N-  1 ) 

= [n  ] Ps 

(N-n^P^0 

; (5.7)* 

5-1 


r 


m 


r (N).  . (N+l), 

Ln  ] = [n  J 


= m0P 


(N+l) 

n+l 

o(N) 


rn(N+1)i 

n+l 


P 

o 


(N+l) 


t>(N) 

“ o 


, (N) 


-m  p 
e o 


SP 


(N) 


n 


p(N) 

n 


3p 


n 

P 


a[<»] 

3P 


3 P 


(N) 


a p 


.i£h 


L? 


(N-i; 


(5.8) 


(5. 9)* 


(5. 10) 


(5.11) 


(5. 12)* 


;5.  13)> 


(5.  14) * 


Note  that  each  of  the  equations  whose  identifying  numbers  are  starred  (*) 
has  a form  identical  to  that  obtained  for  the  finite  population  model 
(parallel  server  inventory  model).  The  fundamental  differences  between 
the  two  models  arise  from  the  fact  that  one  model  assumes  a finite  popu- 
lation, while  the  other  assumes  an  infinite  population.  In  particular,  note 
that  the  P^  for  the  Poisson  Model  corresponding  to  those  states  where 

total  mission  requirements  cannot  be  fulfilled  (n  > N - m +1)  have  de- 

o 

rivatives  and  magnitudes  which  are  very  different  from  those  given  by 
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r~ 


the  Parallel  Server  Model.  Also,  the  Poisson  Model  defines  an  infinite 
number  of  states  associated  with  the  backorder  condition,  while  the 
Parallel  Server  Model  defines  a finite  number  of  states  (m,  in  number). 


SECTION  6 

CONSERVATIVE  SINGLE  SERVER  INVENTORY  MODEL 


The  queueing  model  for  recoverable  items  diagrammed  below  has 
been  investigated,  and  the  steady- state  solution  has  been  obtained. 


The  parameters  for  the  system  diagrammed  above  are  given  as  follows: 

N = total  number  of  items  in  the  system, 
n = number  of  items  in  resupply, 

m = number  ox  items  installed, 

s = number  of  items  in  stock, 

X = failure  rate  for  a single  item  while  installed, 

^■m  = failure  rate  of  installed  items  when  number  of  items  installed 

is  m, 

mo  = minimum  number  of  serviceable  installed  items  required  to 
accomplish  mission  objectives  (m>  m ^ implies  that  service 
performance  level  is  being  met;  m <m  implies  that  the  in- 
stalled serviceable  items  are  working  at  full  capacity  and  are 
not  fully  meeting  mission  requirements), 

Bi  s minimum  number  of  installed  serviceable  items  resulting  in 


The  above  requirement  for  demand  rate  is  equivalent  to 


Xa  - (N  - n)A  , N - m0  + 1 < n S N , 


0 < n < N - m. 


The  state  transition  probabilities  are  given  by 


3Po  ; (t) 


- -*0xpW(t)  + upW(t) 


,00,,. 


3Pv,'(t) 

n 


-tn0XP<N)(t)  + uP^Ct)  - UP^N)Ct)  + maXP^)i(t), 


SP^Ct)' 

u . 


-(N-a;\PW(t)  + uP^(t)  - uP^N)(t)  + (N-n+DXP^Ct), 


3P^N)Ct) 


uP^N)Ct)  + XP^Ct) 


n ■ N . 


In  the  above,  P (t)  is  the  probability  that  n items  are  in  resupply  at 


(6.  6) 


(6.  7) 


(6.8) 


0 < n < H - ; (6.  9) 


N > n > N - m.  + 1 ; (6.  10) 


(6. 11) 


time  t.  Using  standard  techniques,  one  obtains  the  following  steady-state 


solution  (t  = ®). 


I 


(6.  12) 


N-m  m0!  n-(N-o0) 


(6.  13) 


(m0p)  TNI55T  p 


Equation  (6.  13)  can  be  rewritten 


(m0p)N  ma  mo  I “ry  ?oN) 


(6.  14) 


where  x = N - n.  The  probability  that  there  are  zero  units  in  resupply  is 


(6.  15) 


N-m0+l 


N o 


+ (1  - m0p)(m0p)  m0  !G (p) 


where 


G_(p) 


(6.  16) 


Note  that  the  backorder  level  is  zero  for  0<  n<  N - and  is  equal  to 
m x - x = nij  - (N  - n)  for  N-m  + 1 < n<  N.  Thus,  the  expected  backorder 
level  is 


(6.  17) 


so  that 


[N<N)]  - (m1  - x)(m0p)NCni0!mtt  °)P0(N)  ^ (ax  - x)  . (6. 


Solving  for  N = N([N^  ],  p)  obtains 

B 

|-in(m0p)|  |n(P,[N^N)])| 


„ m -l 

“®n  .-X  J— 


— + Znj  (l~m0p)  (m1-x)m() ! m •£— - + 


o ^+21.  (mrx)(a0p)- 

x*o 


*[H«]  C»,pfV‘ 


- m0p)m0  m0!  G_(p) 


(6.  19) 


If  m}  = ij,  then 


[n<n)] 


(oqP)N  mo  ° m0!(l  - m0p)  | ®0Gm<P>  ~ £ Gm0-i(p)} 


1 - (®0P) 


N-m„+i 


N -mo 

+ (1  - map)(m0p)  m9  m0  !Gffl(p) 


» “o  “ “i» 


(6. 20) 


|-dn(«9P)"l  [nCp.Cn^5]  j” 


-6t[Ng  ' ] 


~mo 

mo  mo 


!(1  - m0p)Gn(p){m  - [N*  ']} 


-m.+i 

+ [N^N)](m0p) 


1 *) 
--m0  m0t(l-o0p)GBio_1(p)|  , 


(6. 21) 


®o  " mi  * 


For  a given  fixed  value  of  [N  ' ],  equation  (6.  19)  gives  the  relationship 

of  N and  p.  Thus,  all  values  of  (N,  p)  that  yield  a particular  value  of 
[NjN)]  may  be  obtained  for  comparison  of  inventory  investment  level  to 
resupply  time  level  for  purposes  of  trade-off  analysis  described  earlier 
in  this  document. 

Note  that  when  mQ=  1,  then  Xm  = X for  allm  > 0,  corresponding  to  the 
situation  where  the  failure  rate  is  constant  and  independent  of  the  number 
of  installed  serviceable  units  (except  for  m = 0).  Returning  to  equation 
(6.  15)  and  equation  (6.  18),  one  obtains 


0°  (1  ~ P) 
1 - p1*1 


ma  " 1 


0 < n < N 


(6.  22) 


1 - P 
N+i 


1 - P 


ma  - i 


(6.  23) 


N— uij+i  r ] 

CnBN)]  “ £jpi“  i - p - P)  + i]  j . ®0  - 1 . 

a - pki  - p1^  ) ^ ) 


(6.  24) 


F urther , 


Cm.-D-tno  ] + £t(l-o)  - U 

*«.[■?”»  2 ZZ- 


{ 1 -o“5  [l  ♦ (1-0)  ])]} 


(6.  25) 


«,  • l 


* 


5N  . ££2±ii  «s0  for  [N^]  - con* tan t 


When  N = p = m1=®0  = l,  then 


- \ ; 


3C«i°] 


4 * l 3p 


3p/C^°] 


- 1 . 


Thus,  if  [N  ' ] = constant,  then 

5 


SN  = \4  t ; 

H 


(6.32) 


if  p = constant,  then 


<5Cn4x>3  - 


-t  « ; 


(6.33) 


if  N = constant,  then 


4C4°:  - T* 


4 'W‘U 


(6.34) 


I S[NgN)] ' 

It  is  interesting  to  calculate  [nJN+  - [N^^]  and  — rr in 

B B \ ap  j n 

a fashion  similar  to  that  used  in  the  case  of  the  conservative  parallel 
server  model.  Inspection  of  equation  (7.  12)  and  equation  (7.  13)  obtains 


» ■ - 


J 


,(N+0 


S™.  - KP) 


o < x < n ; 


"N-r+-i 


(6.35) 


and  returning  to  equation  (7.  17)  leads  to 


p(N+0 

" " 0 

ro 


(6. 36) 


Equation  (6.  12)  and  equation  (6.  13)  yield 


/3p00 

1 /3FN-x 


73*T  \3p 


N - x , 1 

P p(») 

0 


0 < x < N ; (6.  37) 


(6. 38) 


Thus, 


WN,]\  [nb,n,i 


(N),  [N  'N'f 

L_  L-»  > - [nIN,i  ♦ > > » • (6.  39) 


an  expression  which  is  identical  in  form  to  that  obtained  for  the  parallel 
server  case. 

Other  formulae  which  apply  are  given  below. 

[n'NU)]  - [n'n>]  , 1[N'N»]  . -p'N+1>  [N<N'l  * 16-40 


SECTION  7 

EXAMPLE  COST  FUNCTION 


The  total  system  cost  is  given  by 


C = C(N,  p)  , 


(7.1) 


where 


N = number  of  inventory  investment  items, 
x T 

p = — = —I - = ratio  of  mean  resupply  time  to  mean  time 

^ X between  demands. 

Total  system  cost  is  expected  to  increase  linearly  with  increasing  N and 
to  increase  with  decreasing  p (decreasing  t ).  For  purposes  of  exposition, 
assume  a simple  cost  function  of  the  form 


where 


Thus , 


c N 

C(N,  p ) = c j N + ~ 


c j = cost  per  investment  item, 

c2  = cost  per  investment  item  per  unit  p. 


inventory  investment  cost  - c N , 


(7.2) 


cost  for  resupply  system  - 


Define 


z(N,  p ) = 


CULP) 


N 1 + 


(7.3) 


P»  = 7“ 


(7.4) 


Thus, 


K 


z(N,  p) 


N 1 + 


(7.5) 


z(N,  P)  is  a unitiess  cost  function  expressed  in  terms  of  number  of  equiva- 
lent investment  items. 


z + z 
1 2 


(7.6) 


Zj  = N = number  of  investment  items  , 

Np0 

z = = cost  of  resupply  system  in  terms 

^ of  equivalent  investment  items 

P0 

Tf  — = 1,  then  for  every  dollar  spent  for  inventory  investment  an  equal 

P 

amount  is  spent  for  the  resupply  system. 


i 
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SECTION  8 

PROBLEM  FORMULATION 


The  mathematical  statement  of  the  problem  takes  on  two  possible 
forms. 


PROBLEM  1 


Minimize 


z = z(N,  p) 


(8.1) 


Subject  to 


[NT,(N,  p)]  = nJ0)  = constant 

N = integer  > 0 
p > 0 

PROBLEM  2 


Minimize 


[nb]  = [Nb(N,P)] 


Subject  to 


z(N,  p)  = zQ  = constant 
N = integer  > 0 
p > 0 


(8.2) 


(8.3) 


(8.4) 
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SECTION  9 

PROBLEM  SOLUTION 


In  this  section  are  presented  the  procedures  to  be  followed  for 
problem  solution.  PROBLEM  1 and  PROBLEM  2 are  restated  below,  and 
associated  optimality  conditions  are  given. 

PROBLEM  1 

Minimize 


z = z(N,  p ) 


(9.  1) 


Subject  to 


P(N.J  \ N)  = P(N),  with  nJ  ^ fixed. 

B B 


(9.2) 


z = N 1 + 


(o  ) 


P(N^  , N)  / 


z = N 


1 + 


P(N) 


Since  N is  an  integer. 


z{N+l,  p(N+l)}  - z{N,  p(N)}  >0 


and 


(9.3) 


(9.4) 


(9.5) 


z{N-l,  p(N-l)}  - z{N,  p(N)J  > 0 , (9.6) 


when  z is  a minimum. 


PROBLEM  2 


Minimize 


[Nb1  = [NB(N,p)]  , 


(9.7) 


Subject  to 


p = p(N,z  ) - p(N)  , with  z0  fixed. 


(9.8) 


P(N)  = NP°  , N < zQ  ; (9.9) 

z -N 
0 

[NB(N,  P)]  = [Nb(n,  p(N)3]  . (9.10) 


If  N is  an  integer,  then 


[N_{N+1,p(N+1)}]  > [N_Cn,P(N)3]  (9.11) 

JD  -D 


and 


[N_  £n-1,  p(N-l)}]  > [N_{N,P(N)}]  , (9.12) 

n>  id 


when  [Ng] 


is  a minimum. 


Both  of  the  above  integer  programming  problems  are  solved  by  a 
simple  search  utilizing  a gradient  technique  (described  in  Section  10)  that 
obtains  solutions  in  two  to  six  interactions.  A similar  technique  is 
utilized  to  obtain  solutions  to  the  equation 


(o ) 

p = p(Ng  , N)  , 


(o  ) 

where  p is  to  be  determined  when  N and  the  backorder  level  NB  are 
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specified.  In  general,  PROBLEM  1 requires  more  computing  time  than 
that  required  for  PROBLEM  2.  The  solution  to  p = p(Ng°*,N)  requires  a 
great  deal  more  time  to  obtain  than  that  required  for  the  solution  of 


z -N 
o 

The  mathematical  programming  problems  described  above  have 
been  thoroughly  investigated  for  the  backorder  function  described  by 
(1)  the  Parallel  Server  Model  and  for  (2)  the  Poisson  Model.  Compari- 
sons and  conclusions  are  given  in  Section  12. 


SECTION  10 

A GRADIENT  TECHNIQUE  TO  SOLVE  A SET  OF  COUPLED  EQUATIONS 


be  n component  column  vectors  and  the  gradient  operator  V be  defined  to 
be  the  row  vector 


dxl  ’ *x2  ' ‘ ' 


It  is  desired  to  find  a solution  to  the  set  of  n simultaneous  equations 


g.(r)  = 0 , i = 1,  2,  3,  . . . , n 


in  n unknowns  represented  by  the  column  vector  x.  The  first  order  Taylor 
expansion  about  a given  point  is 

g.(£  + 5)  = Zg-®  + g.(r  ) 

1 O 1 1 o 

where  Vg  = Vg  (r  ) ; i.e.,  the  first  derivatives  are  evaluated  at  the  point 
i i -o 

r , and  r is  assumed  to  be  an  approximate  solution  to  equations  (10.3). 

“ o — o 

If  r = r + £ is  to  be  a solution  to  equation  (10.  3),  then 
~ ~o 


g.(r_  + e)  = 0,  i = 1,  2, 


V g.  s — ”g.(r)  , i — 1,  2,  . . . , n 

- 1 - 1 ~o 


(10. 6) 


Restating  equation  (10.6)  leads  to 


n 

I (£*,).  «j  • 

J=1  j 


(10.7) 


Now  define  the  square  n x n matrix  G by 


H rH 


I Sx 


'5X2'o 


i<£0> 


2»n(£.7 


(10.8) 


and  the  n component  column  vector  g(rQ)  by 


81(V 

*2<-v 


(10.9) 


g (£  ) 
n o 


In  the  above  it  is  understood  that  all  derivatives  and  all  functions 

V*ih 

g (r  ) are  evaluated  at  the  point  r = r . It  is  observed  that 
si-o  - -o 

[G(r  )]..  = [Vg.(r  )].  and  [g(r  )].  = g.(r  ) • (10.10) 

o 13  1-03  o 1 1-0 


Therefore, 


n 

Y G..(r  )€. 

3 = 1 '3  ‘°  3 


-§i(-ro)  1 


G(r  )«  = -g(r  ) 

— - o — o 


The  above  implies  that 


(10. 11) 


' * -5'  <£  >s(£j  • 


(10.  12) 


-1 

-G  £ 


if  the  matrix  G exists  (i.  e.,  if  det|G|i*0).  The  individual  components 


(10.  13) 


of  £ are  given  by 


n 

X!  Gjkgk  ’ 


j=l,  2,  . . n. 


(10.  14) 


As  long  as  r is  "reasonably  close"  to  the  particular  solution  desired  and 
— o 

no  pathologic  conditions  exist  (forcing  det|G|=0),  then  equation  (10.  13) 
can  be  used  to  successive  arrive  at  values  of  r^  which  satisfy  equation 
(10.3)  as  precisely  as  desired. 

The  procedure  simply  stated  is  as  follows: 

(i)  select  an  initial  trial  solution  value  £ = r^,  where  i = the  iteration 
number  (i  = 0 initially); 
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(ii)  apply  equation  (10.  13)  to  determine  the  correction  vector  e.; 

— 1 

(iii)  check  to  see  if  e.  satisfies  I e.  |<  e , where  the  (e  ).  = 

1 ‘ — i -max  -max  j 

(e.)  >0  represent  upper  limit  values  on  the  errors  to  be 

j max 

tolerated  in  the  solution  values  for  each  of  the  variables  x.; 

J 

(iv)  (a)  if  I e.  I > € , then  set  r.  , = r.  + e.  (i.  e.  , calculate  a new 

‘“l  -max  ~i+l  -i  -i 

trial  solution  to  the  system  of  equations  g(r)  = 0)  and  return  to 

step  (ii)  above  with  i = i + 1;  (b)  if  I e.  |<  e , then  stop  and 

-i  —max 

calculate  r.  , = r.  + e.,  where  r.  , represents  the  desired  solu- 
-l+l  -l  -i  -l+l 

tion  to  g(r)  = 0. 

As  an  example,  consider  the  problem 

2 2 

gx(x,y)  = x - y - 4 = 0 

g2(x,  y)  = x + y - 1 = 0 

The  solution  to  the  above  problem  is  readily  obtained  by  conventional  tech- 
niques: 

x + y = 1 ; 

x - y = 4 ; 

2x  = 5 

x = 5/2  and  y = -3/2 

To  illustrate  the  technique  developed  above  proceed  as  follows: 


4g 


1 


(2x,  -2y) 


2x  -2y 1 


= (1,  1) 


G = 


det  G 


2(x  + y) 


2(x  + y) 


2y' 


2x> 


g = 


, 2 2 
x - y - 4 


x + y - 
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Choose  r 
“o 


as  the  initial  trial  solution 


y W 


Iteration  No.  1 proceeds  as  follows: 


G-1(r  ) 


tu  | oj  | m 


Note  that  it  is  not  necessary  that  the  first  approximation  be 
terribly  good;  what  is  required  is  that  the  chosen  be  "closer"  to  the 
solution  of  interest  than  to  any  other  possible  solution  (thus  preventing 
the  possibility  of  moving  along  in  the  wrong  direction)  and  that  in  succes- 
sive calculations  no  pathological  conditions  causing  det|G|  to  be  zero 
occur  [such  as  encountering  a saddle  point  or  a relative  minimum  or 
maximum  for  any  of  the  g^(x)] . If  such  a condition  does  arise,  however, 
it  is  still  possible  to  apply  the  technique  by  "skipping  over"  the  critical 
point  encountered  [Note  that  det|G|  = 0 for  any  x which  is  a critical  point 
for  any  of  the  g (x)  in  g(x).  ].  A "good"  value  for  tq  eliminates  such  a 
possibility  from  occurring,  and  in  general  no  such  problem  occurs  in 
practical  applications. 

The  one  dimensional  case  of  the  solution  technique  is  a special  case 
of  the  development  presented  above.  To  obtain  a solution  to  the  equation 

g(x)  = 0 ( 10.  15) 

note  that 

Jdg(x)\ 

g(x  + e)  = — c + g(x  ) (10.16) 

8 O \ dx  / o ° 

and 

-1 

' * -s'V  (If)  0 ,10-17) 

if  x = x + e is  to  be  a solution  to  equation  (10.  15).  The  procedure  to 
o 

determine  the  solution  is  precisely  the  same  as  that  given  in  steps  (i) 
through  (iv)  above. 
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\ 


An  example  Illustration  of  the  one  dimensional  case  is  given  below. 
In  order  to  solve 

x + 2nx  = 0 

for  x,  one  proceeds  as  follows: 

g(x)  = x + 4nx  ; 


dx 


i ♦ i 

x 


x(x  + Inx) 
1 + x 


Try  x = 1 as  a trial  solution.  Then 


e = -0.500. 

o 


Iteration  No.  1 proceeds  as  follows: 


x,  = x + e =0.5  ; 

1 o o 


k" 

i+i 


T 1. 5 

Iteration  No.  2 proceeds  as  follows: 


0.  5(0,5  - 0.  693147).  = +0.0643824 


X2  = Xl  + €1 


0.  5 + 0.0643824  = 0.  5643824  ; 


[0.5643824)  (0.5643824  - 0.5720233) 
1. 5643824 


e = +0.00275659 

u 


Iteration  No.  3 proceeds  as  follows: 


x = x +c  ^ 0.5671389835  ; 

D fa  b 


■ (0.5671389835)  (0.5671389835  - 0.5671508844) 
1.5671390 


+0. 0000043069 
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Iteration  No.  4 proceeds  as  follows: 


x = x + c,  = 0.5671432904 

4 3 3 


-(0.5671432904)  (0.  5671432904  - 0.5671432904) 
®4  = 1.5671432904 


x = 0.  5671432904  to  nine  significant  figures 


If  a solution  value  to  the  equation  g(x)  = 0 is  desired  when  x is  con- 
strained to  be  an  integer,  then  the  above  procedure  is  modified  as  follows. 


g(x  + e)  = g(x  ) + Ag(x  )e  = 0 

O o o 


(10. 18) 


Ag(x  ) = g(x  + i; 

o o 


g(xo) 


(10.  19) 


[- 


g(xo>1 

^xo) 


(10. 20) 


where  the  subscript  T means  set  e = integer  = truncated  value  of  the  numeri- 
cal quantity  enclosed  in  the  brackets.  The  iterative  procedure  is  then  that 
given  by  steps  (i)  through  (iv)  given  earlier  in  this  section  with  equation 
(10.20)  being  substituted  for  equation  (10.  17)  and  the  initial  value  for  xq 
being  set  equal  to  an  integer.  If  no  solution  to  g(x)  = 0 exists  (this  is  possi- 
ble since  x = integer),  then  the  above  procedure  may  be  utilized  in  finding 
the  value  of  x .>  0 that  minimizes  j g(x)  | and  the  value  of  x < 0 (if  allowed) 
that  minimizes  j g(x)  | . 

In  the  solution  procedure  for  PROBLEM  1,  the  gradient  technique  is 
applied  to  determine  p in  equation  (9.  2)  and  is  used  to  aid  in  the  search  for 
the  optimal  value  of  N which  satisfies  the  conditions  expressed  in  relations 
(9.5)  and  (9.6).  The  search  for  (N,  p)  is  conducted  along  a curve  of  constant 
[N  ].  In  the  solution  procedure  for  PROBLEM  2,  the  gradient  technique  is 
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used  to  aid  in  the  search  for  the  optimal  value  of  N which  satisfies  the 
conditions  expressed  in  relations  (9.  11)  and  (9.  12).  Here,  the  search 
for  (N,  p)  is  conducted  along  a curve  of  constant  z. 

The  solution  procedure  developed  here  for  PROBLEM  1 and  PROB 
LEM  2 is  easily  adapted  to  any  cost  function  of  the  form  C = C(N,  p)  and 
for  any  backorder  function  of  the  form  [N  ] = [N  (N,  p)].  In  this  study, 
PROBLEM  1 and  PROBLEM  2 have  been  investigated  with  the  backorder 
functions  described  by  the  Poisson  Model  (Section  5)  and  the  Parallel 
Server  Model  (Section  4)  and  with  z = z(N,  p)  as  described  in  Section  7. 
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SECTION  11 

SENSITIVITY  ANALYSIS 


It  is  of  interest  to  determine  the  range  of  p for  which  the  optimum 
value  of  N remains  unchanged.  For  PROBLEM  1,  application  of  the  opti- 
mality condition  and  the  corresponding  constraint  yields  values  for  (p„  )max 

and  (p  ) such  that  when 
o mm 


(p  ) . < P < (P  ) 

o min  o o max 


(11.  1) 


the  current  optimal  value  obtained  for  the  (N,  p)  ordered  pair  remains  un- 
changed; i.  e.  , the  optimal  solution  to  PROBLEM  1 is  independent  of  the 
choice  of  p within  the  interval  (P0)min^  PQ  ^ ^o^max’  Similarly,  appli- 
cation of  the  optimality  condition  for  PROBLEM  2 and  the  corresponding 

constraint  yields  values  for  (p  ) and  (p  ) such  that  when 
1 o max  0 min 


(P0> 


mm 


^ P. 


<Vma* 


(11.2) 


the  current  optimal  value  of  N remains  unchanged;  i.e.  , the  optimal  solu- 
tion value  for  N for  PROBLEM  2 is  independent  of  the  choice  of  pQ  within 

the  interval  (pj  ^ p.  ^(P„)  . • However,  in  the  case  for  PROBLEM  2 

0 max  0 0 mm 

if  p is  allowed  to  vary  and  the  optimum  value  of  N remains  fixed,  then 
o 

the  optimum  value  of  p must  change  to  keep  z = constant. 

Let  (N,  p)  be  an  optimal  solution  to  PROBLEM  1.  Define  Ap:  ^ 0 and 
Ap2  > 0 by 

[Nb(Ni  + 1»  P , + Ap ! )]  = N^°}  , (11.3) 

[N  (N  - 1,  P - Ap  )]  = N*0)  . (11.4) 
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Condition  (9.5)  in  Section  9 becomes 


(Nx  + 1)  1 + 


5 N. (1  + p./p.>  (U-5) 


and  yields 


1 + APj  /px 


p./p!  £ — at 


N,  — - 1 


Condition  (9.6)  in  Section  9 becomes 


(11.6) 


,N>  - 11  1 * P - N>  (1  + po /pi  > 


(11.7) 


and  yields 


1 - Ap  /p 
2 1 


(11.8) 


Thus,  any  P0  satisfying  the  condition 


1 

P 


< p / p < 

— K0  ' H1  — 


N - 1 

1 Pi 


1 + Ap  j/pi 


N — - - 1 
1 P! 


(11.9) 


will  yield  an  optimal  solution  (N^  , p ) for  PROBLEM  1. 

Let  (Nj,  Pj)  now  be  an  optimum  solution  to  PROBLEM  2.  Define 

Ap  >0  and  Ap  > 0 by 
l 2 1 


[N„(N  + 1,  p + Ap  )]  = [N  (N  ,p  )]  , 

B : l 


(11. 10) 
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Conditions  (11.  9)  and  (11.  16)  can  be  utilized  to  study  the  sensitivity 
of  the  optimal  fractional  investment  level  to  the  relative  unit  cost  associ- 
ated with  the  resupply  system.  Results  of  comparing  the  Poisson  Model 
to  the  Parallel  Server  Model  are  given  in  Section  12. 
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SECTION  12 

RESULTS  AND  CONCLUSIONS 


Computer  programs  have  been  written  and  executed  for 

1)  The  Conservative  Parallel  Server  Inventory  Model 
described  in  Section  4, 

2)  The  Poisson  Model  described  in  Section  5, 

3)  The  Conservative  Single  Server  Inventory  Model  des- 
cribed in  Section  6. 

Programs  have  also  been  written  and  executed  for  solving  PROBLEM  1 
and  PROBLEM  2 described  in  Sections  8 and  9 utilizing  the  gradient 
technique  described  in  Section  10  and  procedures  outlined  in  Section  9. 
Included  in  these  programs  are  the  sensitivity  analyses  given  in  Section  11. 
Program  listings  are  found  in  Appendix  A. 

Reference  material  developed  during  the  course  of  this  project  is 
given  in  the  "Reference  List  and  Bibliography"  and  in  "Documents 
Internal  to  AFLC". 

Various  inputs  were  assumed  to  execute  the  computer  programs  to 
perform  a detailed  analysis  and  comparison  of  the  behavior  of  the  parallel 
server  model  (finite  population)  and  the  Poisson  model  (infinite  popula- 
tion). Below  are  listed  the  major  results. 

• The  Poisson  demand  process  overestimates  the  backorder  level 
for  a given  (N,  p ) specification. 

• In.  general,  the  Poisson  Model  obtains  a lower  value  for  the 
optimum  value  of  inventory  investment  than  does  the  finite  popu- 
lation model;  as  the  system  requirement  mQ  increases,  the 
difference  between  the  optimum  levels  of  inventory  investment 
given  by  the  two  models  becomes  increasingly  larger.  Thus, 
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use  of  the  Poisson  Model  underestimates  optimum  inventory 
investment  level  and  results  in  a higher  backorder  level  and  a 
higher  stockout  probability  than  does  the  finite  population  level. 

• If  the  relative  unit  cost  (pQ)  for  resupply  increases,  then  the 
difference  between  the  optimum  levels  of  inventory  investment 
given  by  the  two  models  becomes  increasingly  larger,  approach- 
ing levels  of  15-50%. 

• In  general,  the  finite  population  model  given  an  optimum  level 
for  inventory  investment  which  is  less  sensitive  to  a change  in 
unit  resupply  cost  than  that  suggested  by  the  Poisson  Model; 

i.  e.  , the  range  for  pQ  where  the  optimal  solution  remains  invari- 
ent  is  larger  for  the  finite  population  model  than  for  the  infinite 
population  model. 

• As  the  system  requirement  becomes  larger,  the  optimum 
level  for  inventory  investment  becomes  more  sensitive  to  a 
change  in  unit  resupply  cost. 

Some  numerical  outputs  are  tabulated  in  Appendix  B. 

The  above  results  conclude  that  if  the  relative  cost  for  resupply 
(manpower,  equipment,  administration,  transportation,  etc.  ) compared 
to  inventory  investment  increases,  then  it  becomes  increasingly  impor- 
tant that  demand  and  resupply  processes  be  acc  irately  described  to 
obtain  cost  effective  allocations  of  expenditures  between  inventory 
investment  levels  (capital  intensive)  and  resupply  time  levels  (labor 
intensive).  In  addition,  if  the  variance  (uncertainty)  of  costs  increases, 
then  it  is  also  important  that  correct  descriptions  of  the  demand  and 
resupply  processes  be  utilized  in  the  decision  making  process.  It  is 
noted  that  inventory  investment  costs  for  recoverable  items  are  well- 
known  in  comparison  to  projected  costs  for  resupply  during  the  life 
cycle  of  a typical  weapons  system.  Thus,  the  fact  that  the  finite  popu- 
lation model  tends  to  allocate  a greater  fraction  of  total  expenditures  to 
inventory  investment  than  does  the  infinite  population  model  makes  the 
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finite  population  model  even  more  attractive  for  adoption  in  the  decision 
making  process. 

Non- optimal  allocations  that  result  when  incorrect  demand  and 
resupply  processes  are  assumed  suggest  the  following  areas  for  future 
investigation. 

• Study  resupply  processes  which  allow  for  non- zero  queue  lengths 
(finite  number  of  parallel  servers). 

• Generalize  the  trade-off  analysis  procedure  developed  in  this 
report  (resupply  time  level  versus  inventory  investment  level 
to  minimize  total  system  cost)  to  include  interactions  due  to 
movement  of  many  investment  items  in  the  transportation  pipe- 
lines among  user  bases  and  supply  depots. 

• Further  investigate  demand  and  resupply  modeling  for  recoverable 
items  for  purposes  of  determining  transportation  flows  and  optimal 
allocating  of  resources  to  inventory  investment  level  and  resupply 
time  level. 

• Investigate  and  compare  the  utility  of  the  finite  population  model 
developed  in  this  report  with  that  of  the  Poisson  model  currently 
adopted  by  Air  Force  management  [88,  89,  104]  when  incorporated 
into  the  METRIC  and  MOD- METRIC  recoverable  item  control 
and  allocation  models. 


• Investigate  and  compare  the  finite  population  model  with  the 
Poisson  model  in  determining  impacts  on  output  of  the  Compre- 
hensive Engine  Management  System  (CEMS)  being  developed  by 
the  U.  S.  Air  Force  [D-ll]. 

The  above  areas  of  investigation  and  those  described  in  this  report  should 
contribute  to  a better  understanding  of  certain  problems  of  immediate 
concern  to  the  U.  S.  Air  Force.  The  questions  dealt  with  here  are  very 
important  for  considerations  that  must  be  made  when  an  integrated  design 
effort  is  made  to  develop  any  proposed  weapon  system. 


I 


The  work  reported  in  this  document  was  presented  at  two  technical 
conferences  which  took  place  during  the  months  of  January  and  February 
of  1978.  Particulars  are  given  below. 

1)  "Analysis  of  Inventory /Re supply  Models  in  Optimal  Inventory 
Investment  Allocation,  " ORSA/TIMS  Special  Interest  Conference 
on  Capital  Budgeting,  " January  30,  31,  and  February  1,  1978, 
Phoenix,  Arizona. 

2)  "Analysis  of  Inventory /Re supply  Models  in  Optimal  Materiel 
Distribution  System  Design,  " ORSA/ TIMS/ AIIE  Special  Interest 
Conference  of  Distribution  Models,  February  20-22,  1978. 

Hilton  Head  Island,  South  Carolina. 

The  above  conferences  were  sponsored  by  the  Operations  Research 
Society  of  America  (ORSA),  The  Institute  of  Management  Sciences  (TIMS), 
and  the  American  Institute  of  Industrial  Engineers  (AIIE). 
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APPENDIX  A 
PROGRAM  LISTINGS 


CONSERVATIVE  PARALLEL  SERVER  INVENTORY  MODEL 


(FINITE  POPULATION)  PROGRAM  LISTING  A- 2 

POISSON  MODEL  (INFINITE  POPULATION)  PROGRAM 

LISTING  A- 17 

CONSERVATIVE  SINGLE  SERVER  INVENTORY  MODEL 

(FINITE  POPULATION)  PROGRAM  LISTING  A- 27 

PROBLEM  1 PROGRAM  LISTING  A-45 

PROBLEM  2 PROGRAM  LISTING  A- 5 2 


CONSERVATIVE  PARALLEL  SERVER 
INVENTORY  MODEL 
(FINITE  POPULATION) 
PROGRAM  LISTING 


/LOGON  “NC" 


/DO  RUN«WATFIV 

» JOB  <L8NTJQQ> KOVACS 

DOUBLE  PRECI  SION  PC  1 01  >«  PAC  101  >«DBB2C  1 01  )#NB#NB1.» PO# PO  UZ»U1 

1MUM0»ANS19»  PERU  PER 3*  VAR,  VARA#  NB  SO,  VARN# RES  P SO#  VARS#  SM#MN#  VABM#i) 

2SSQ#MNSQ#T  

Z-0, 

PRINTS 

5 FORMATC  * 1 S 5X#  *»»»KOVACS  FININITE  SERVER  INVENTORY  HO 

BEADC  5#»>  NT0T#R0#M0#M  1 

PRINT  1 0#NT0T#B0#MQ#M1 

1J> FORMATC  * ■ NTQT»,»D20«  14#  * RQ«»*#D20»  14#  D20.  14#*  . Ml««# D20, 14) 

ITOT-NTOT 


PRINT 150 


160  FORMA TC • OTHE  TOTAL  OF  ALL  THE  P VALUES  IS  *# D20#  i 4)  _ 

- 

PRINT 1 60* T 

CALL  S1XTENCNT0T#M1#F#NB) 

PRINT 140# MB  

I6-NT0T-M1+1 j 

IFCI6.LE#0)Z-P0  

I FC  I 6#LE«  0 > 2 6-1 

DO  300  1 7- 1 6#  I TOT | 

3 00 Z-Z+PC  175 

C Z 1J_THE  TOTAL  OF  ALL  THE  PROBABILITIES  FROM  I-NTOT-M  1*  1 TO  I-MTDT#. 

, 

PRINT  301#  Z 

3 01 FORMA Tt  ■ OTHE  PROBABILI TY  FOR  STOCKOUT  I S _•#  D20#14> 

2X-IY-0  , j 

- - - ; 

CALL  N1NE<NT0T*1#M0#PA#P01#R0#IX#IY>  j 

_ C ALLTENC NTO  T+  1#M0#  PA#P01#R0#  IX#  2 Y> 

179-IT0T+1  

PRINT220# IX# I Y 

2 20_  FORMATC  * OTHERE  WERE  *#I3#  * • UNDERFLOWS  AND  *#1  3#  OVERFLO  WS*) 


DO  102  1-1# 279 

1 02 T-T*PA<  1 ) 

PRINT  1 60# T ....  

CALL  SIXTENCNTOT*  1#M1#PA#NB1  > _ 

PR2NT1 40# NB 1 

PRINT  1 42#  NB-NB 1 

1 42  FORMATC*  . DELTA  NB-*#D20#  14) 

PR2NT400# I TOT# I TOT* 1 # I TOT# ♦< NB-NB i 1 ) /NB 
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4 00  FORMAT!  * ONB!  • # I 3#  * > - NB!*#2  3#*>  / NB!*#I3#*>  • *#D20.14) 

CALL  SUB2KNT0T# P#  DERI  > 

PRINT3 1 0#  DERI 

3 10  FORMAT!  • OTHE  EXPECTED  NUMBER  OF  UNITS  IN  RESUPPLY  IS  SD20#14> 

PRINT32  0#  DER 1 /<  0.0- RO > 

3 20  FORMAT!  • OTHE  RESULT  OF  FORMULA  21  IS  *#D24.14)  - ' “ 

- - L.  ~ 

CALL  SUB20!  NTOT# RO* DERI#  DER2) 

PR2NT333 

PRINT330#  0#  DERI /! 0.* RO  > 

PRINT330#! 2# DER2! I ># I- 1# I TOT) 

3 33 FORMAT!  • OANSWERS  FOR  FORMULA  20*  • ) 

3 30  FORMAT!  * ■ P!  *#  1 3#  * ) GIVES  *#D24.14) 

CALL_VAR 1 ! PO#  M 1 # NTOT#  VAR# P#NB> 

CALL  UAR1!P01#M1#NT0T+ 1# VARA#  PA#NB1 ) 

PRINT350# 2 TOT# VAR 

3 50  FORMAT!  * ■ THE  VARIANCE  OF  NB  FOR  NTOT- •#  1 3#  * 2S_*#D20.  143 

PR2NT350#  2 TOT*  1#  VARA 

DER 3-0.0 

CALL  VARK  P0#M1#NT0T#NBSQ#  P#  DER3)  

CALL  SUB22!NT0T#M1#NB#  DERI#  DER3# R0#NBSQ) 

PRINT360# DER3  

3 60  FORMAT!  ••THE  ANSWER  TO  FORMULA  22  IS  •#  D20.  1 4) 

PHI  NT 352#  NBSQ 

CALL  VAR2!  DER  1#P#VARN#P0# NTOT) 

3 52  FORMAT!*  OTHE  EXPECTED  VALUE  OF  THE  SQUARE  OF  THE  BACKORDER  LBVELI  I 
1 S! *#  D20. 1 4) 

PRINT41 0# VARN 

I 


A-  5 


4 10  FORMATC  • OTHE  VARIANCE  OF  THE  EXPECTED  NUMBER  OF  UNITS  IN  RESUPH.Y 


1 1 St  •#  D20«  1 4> 

DER3"0*0 

CALL  VAR2C  DER3j  PjRESPSQ*  PQjNTOT?  

PR  2NT420#  RE  SPS  Q 

4 20 FORMATC  * OTHE  EXPECTED  VALUE  OF  THE  SQUARE  OF  THE  NUMBER  OF  UNITSI 

IN  RESUPPLY  I St  D20» 1 4?  

CALL  SNAVECNTOT#MU  P#  P0»  SN> 

PRINT430# SN 

4 30  FORMATC  ■OTHE  EXPECTED  STOCK  LEVEL  IS  •»D£0*14> 

CALL  VAR3CNT0T»MUP»P0»VARS»  SN > 

PRINT440#  VARS 

440 FORMATC  * OTHE  VARIANCE  OF  THE  STOCK  LEVEL  IS  •#D20,14> 

. u-o«o 

CALL  VAR3CNT0T* M U P#PO*  S SQl*U> 

4 25 FORMATC  * OTHE  EXPECTED  VALUE  OF  THE  SQUARE  OF  THE_STOCX  LEVEL_ISt| 

1D20«14) 

PRINT 42 5# SSQ 

CALL  MNA VEC  MTO  T#M  1 # P#  PO#  MN  J 

PRINT 43 5# MN 

4 35  FORMATC  • OTHE  EXPECTED  NUMBER  IN  SERVICE  I S •*D20.  1 4> 

CALL  VAR4CNT0T*M1 »MN» P» VARMj PO)  

PRINT444*  VARM __ 

4 44  FORMATC  * OTHE  VARIANCE  OF  THE  EXPECTED  NUMBER  IN  SERVICE  IS  «#D£  0* 

14>  

CALL  VAR4CNT0T#M1#U#P#MNSQ*P0> 

PRINT446#MNSQ 


A-  6 


446  FORMATC  • OTHE  EXPECTED  VALUE  OF  THE  SQUARE  OF  THE  NUMBER  IN  SERVIB 


1 ISt  *# D20* 1 4) 


STOP 


SUBROUTINE  NINECNTOT# MO#  P#  PO# RO# IX# 1Y) 

C FINDS  THE  PROBABILITIES  FOR  J-0  TO  J-NTOT-MO  (FINDS  P(J>>« 


C USES  NATURAL  LOBS  AND  EXPONENTS  TO  ALLOW  FOR  LARGER  SY STEMS 

DOUBLE  PRECISI  ON  P(  101  ># NTOT#  PO#  MO#  RO#  A#  B#  C#  D#  E#  F#  b#K#  Z#  Y 

DOUBLE  PRECISION  DL06#  DEXP 


ITEST-NTOT-MO+ 1 

DO  10  Ul- 1# I TEST 


Y-Z»0#0 


DO  20  N 1-1# I TEST 
NRESUP-N1-1 


A— < NRESUP- EL  0<3<  M0»  RO  > 


C-D-0.0 


IFC J*LE* 1 2G0T0  35 


^ r\  •“ 
yj  7 


DO  30  X- 2#J 


3 0 C-C*DLOG(B> 


35  IFCNRESUP.LE#  UQOTO  45 


DO  40  K-2# NRESUP 


4 0 E—  D*  ELOG<  B> 

4 5 E-A+C-D 


CALL  FLOVC  E# J# IX# I Y) 
Z-Z+DEXPCE) 


A- 7 


20 


CONTINUE 


50 


CONTINUE 


IFCJ.EQ«Q>PO-UO/<Y+Z  ) 

IF<J.EQ  •OJPRI NT503*  J#  PQ 

IF<  J.EQ.OQOTO  1 0 

P<J5-l./<  Y+Z  ) 

PRINT503* J#  P<  J> 

5 03  FORWATC  » P<  ■ *13#  »>  - ■*DgQ#M> 

1 0 CONTINUE 

RETURN 

END 

SUBROUTINE  TEN<NT0T*M0#  P»  PO*RO#  IX#  IY> 

C FINDS  THE  PROBABILITIES  FOR  J-NTOT-MO+i  TO  J-NTOT  CUSES  NATURAL  ELOQST 

0 ALLOW  FOR  LARGER  NTOT>._ 

DOUBLE  PRECI  SION  P<  101)j SITGT#  P0*H0# RO# A# B# C*  D#  E*  F#  Q»H# Z#  Y 

DOUBLE  PRECISION  ELOQ*  DEXP  

INITAL-NTOT-MO-M  


1 05  N-NTOT-J 


X F<N*LE# 1 5 GOTO  115 

DO  IIP  K-2#N 

B-K 

1 10  H-H*ELOGCB> 


1 

1 15 

IFCNRESUP*LE#  13Q0T0  35  . 

DO  30  K-2#NRESUP  __  „ 

B-K 

30 

A-A+EL0Ci<  B)  _ . . 

35 

N-MO 

I F<N*LE« 1 3GOTO  45 

DO  40  K-2*N 

B-K  „ . 

40 

C-C+CL0G<B> 

4 5 D«<NRESUP-J3*EL0G<  R0> 

E-<  NRESUP-<  N TOT-MO 3 >* EL06< MO  > 

F-P»E»Q»H-A-C 

CALL  FLOW<  F#  J#  I X # I Y > 

Y-Y*DEXP<F>  

2 0 CONTINUE 

DO  50  NRESUP-INITAL#1T0T 

A-C-D-E"F-0#0 

IF<J#LE»1)00T065 

DO  60  K*2»J 

B-K  

60  A-AfELOQ<B>  


A-  10 


6 5 X-NTOT-J 

IF*  2 «LE«  1 >QOTO  75 

DO  70  K»g<I 

B-K 

PPELOQCB? 

i-NTOT-MBESUP 

IF<1<LE<OGOTQ  85 

DO  8 O K- 2# I 

B-K 

E-E+JDLOGCB? 

if<nresup<le<  ugoto  95 

DO  90  K-2<MBESUP 

B-K 

F-F+ELOGC  B) 

D-<NRESUP-J3»EL0Q< BO  3 

B-D*A-»C-E-F 

CALL  FLOVC  B<  J<  IX*  I Y> 

Z-Z+DEXPCB5 

5 0 CONTINUE  

P<  J5-  U0/<  Y+Z  > 

PRIMT506<J<P<  J> 

5 06 FORMATC  • . P<  •<  1 3<  ' ) - •<  D20< 1 4> 

JO  COMTIMUE 

RETURN 

END 


70 

75 


A-  12 


C _ FINDS  THE  RESUL TS  FOR  FORMULA  20*  . 

DOUBLE  PRECISION  NTOT.RO# DER1» DER2C 101> 

ITOT-NTOT  7 ; 

DO  10  NRESUP*1#IT0T  __  

10 DER2C NRESUP>*NRESUP/RO«  DERI /RO  

P ODER* PER  1 /<  0.  O-ROJ 

RETURN 

END 

SUBROUTINE  SUB21CN TO T# P»DER  1 > 

C FINDS  THE^ DEXPECTED  NUMBER  OF  UNI TS  IN  RESUPPLY. 

DOUBLE  PRECI SION  NTOT# P< 101 )• DERI  _ 

DERl-0.0  __  

ITOT-NTOT 

DO  10  NRESUP* 1. I TOT 


A-  12 


H 


10 DER  1»NRESU  P*  P<  NRE  SU  P > ♦ DEB  1 

RETURN  

END^ 

SUBROUT  I N E_SUB  2 SCNTOT#  M l # N B#  DER 1 # DER2#  RO# NBSQ ) 

C FINDS  THE  RESULTS  FOR  FORMULA  22.  . 

DOUBLE  PRECISION  NT0T#M1#NB# DER 1#DER2* RO#NBSQ_ 

DER2»NB/R0»<  (NT0T-M1 ?-DERl+NBSQ/NE> 

RETURN 

ram 

SUBROUTINE  VAB3<NTOT# M 1 # P# PO# VAR#  SN ) 

DOUBLE  PRECISION  NTQT*Ml#P<  1 0 1 >#  PO#  VAR#  SN 

I TES  T»N TOT-M 1 

VAR«P0»<NT0T-M1-SN?»»2  

DO  10  NRESUP»  1#  I TEST  

10 yAR«VAR*P<NRESUP>*<NT0T-Ml-NRESUP-SN5**2 

RETURN 

END 

SUBROUTINE  SNAVE<NTOT»Ml#  P#  PO#  SN>  

DOUBLE  PRECISION  NT0T#H1 # PC  1 01 ># PO# SN 

ITEST-NT0T-M1_ . 

SN-P0*<NT0T-M1> 

DO  10  NRESUP" 1 # I TEST 

SN-SN4PCNRESUP>_*<NT0T-M1-NRESUP>  

1 0 CONTINUE  _____ 

RETURN  ...  

END 


A-  1 3 

L i 


J 


SUBROUTINE  MNAVECNTOT#Ml#P#PO#HN>  . . - 

_ DOUBLE  PRECI SION  NTOT»Ml»PC  101>»PQ»»W 

INITAL-NT0T-M1__ 

MN-0 

I F<  INI  TAL.EQ#  Q?MN-PO»C  NTOT> 

ITOT-NTOT 

IFC  INITAL*  EQ#0JINI  TAL-1 

DO  10  NRESUP- INI TAL# I TOT 

— ■ 

10 MN-MN*<NTOT-NRESUP>* PC NRESUP?  

MN-MN*Ml*PO 

I TEST-INI TAL-1  _ 

_ — ^ 

DO  20  NRESUP- 1*1  TEST 

2 0 MN-MN^M1» PC NRESUP?  j 

RETURN 

END 

SUBROUTINE  VAR4C  N TOT#  M 1 # MN#  P#  VAR#  P0> ___ 

DOUBLE  PREC 1 SI ?N_NTO  T# Ml # MN#P<  101?# VAR# PO 

INITAL-NT0T-M1 

I TEST-INI TAL* 1 

VAR-PO  

DO  20  NRESUP- 1#  I TEST 

2 0 VAR- VAR+PC  NRESUP? 

VAR-  VAR»<  M 1 -MN  >»»2 

IFC  INITAL*EQ*  0>VAR-VAR*PO*CNTOT-MN?**2 

IFC  INITAL#EQ*0?INJ TAL-1  _ 

ITOT-NTOT  

DO  10  NRESUP— INI TAL# I TOT 


A-  14 


» 


I Ft  INI TAL* EQ*  0 J VAR«VAR> PO*<M 1 -NTOT-NB3**2 


IFCJNITAL.EJl.  OJlNITAL- 1 — i 

DO  10  NRESUP-INITAL#ITOT  f 



1J3 UAR"VAR+P<aiRESUP3*<H  1 -NTOT+NRESUP-NB 3**2 

RETURN  


END 

SUBROUTINE  VAR2C  PERU  P*  VAR*  PO»NTOT> 

DOUBLE  PRECISION  P<  101  PERU  VARjPO+NTOT 

ITOT-NTOT 

VAR- PO*DER 1**2 <; 

DO  10  NRESUP-UITOT  i 

1 0 VAR-VAR*P<NRESUP3*<NRESUP7DER1J**2 

— — - • ■ 

RETURN 

END ; 

SENTRY — 

20  *5  15  IT  _ __  ! 

/LOGOFF  I 


/LOGON 


"NC" 


/DO  RUN»WATFIV  

S JOB  C LSNTJQQ > KOVACS  

DOUBLE  PRECISION  PC  101>»PAC  101 )#  DER2C  1 01  )#NB#NB1# FQ#PQ1#  STOCKISTS 

IK  1 # N TO  T#RO# M 1 # M 0#AN  S 1 9#  DER 1 # DER  3#  VAR#  VARA#  MBS  Q#  VA RN#  RES P SO#  VARS#  31 

2#MN#  VARM#  U#  SSQ#MNSQ#  T ■ 

READ#NTQT#BO#  MO#  Ml 

ITOT-NTOT  

CALL  PVALC  P»MTQT#MO#  RO#  PO) 

CALL  NBSTK<NB<STOCK< N TOT# M 1#M 0#  RO#  PO) 

CALL  PVALC  PA#NT0T+1#M0# RO#  PO  1 5 

t 

CALL  MBSTKC  NB 1#  STOCK  1#NT0T»  1#M1#M0#  RO#  P01> 

___  PRINT 142#NB-N  B 1 

CALL  N2NTENCNT0T#M1#NB1#NB#  P01#  PO#  AN  SI  9#  MO#  FA#RO) 

142  FORMATC  • DEL  TA  NB-  • # D2  0 • 1 4 > 

PRINT4 00#  1 TOT#  I TOT*  1#  1 TOT#  ANSI 9 . 

4 00  FORMATC  'CNBC*#  I 3#*)  - NBC  •#23#*>  / NBC«#23#  •>  - *#D20#14> 

PRINT304#*CNE-NB1)/NB  _ 

3 04_  FORMATC'^THELEFT  SI  DE  OF  THE  EQUATION  GIVES ...  *#  D20.  14> 

CALL  SUBS  1 C NTOT#  P#  DER  1 > 

PRINT3i0#DERl  

3 10 FORMATC  «OTHE  EXPECTED  NUMBER  OF  UNITS  IN  RESUPPLY  I S«#D2G«14) 

PRINT320#  DERI /C-RO)  _ 

3 20  FORMATC  • OTHE  RESULT  OF  FORMULA  2 1 I S •#  D24*  1 4> 

CALL  SUB20CNT0T#  RO#  DERI#  DER2>  — - 

PRINT333 


A-  18 


PRINT330#  0#  DERI /C -RO) 


3 33  FORMATC • CANSWERS  FOR  FORMULA  20*  •> 


3 30  FORMATC  * ■ PC  •#  I 3#  • > 02  VES  •#  D24.  1 4 >_ 

PR1NT3 3 0#  CJ  6#  DER2C  J63#J6«1#IT0T>  

CALL  VAR  1C  PC* Ml»NTOT#VAR#  P, NB  > 

CALL  VAB1<P01#M1#NT0T*1#VARA#PA#NB1>  

PRINT350#  2 TOT* 1#  VARA 

f 

350  FORMATC*  THE  VARIANCE  OF  NB  FOR  NTOT»*# 1 3#  * • I S »#D20«14> 

PRI NT35  0#  I TOT+1#  VARA  

DER3-Q.  0 

CALL  VAR1CP0#M1#NT0T#NBSQ# P# DER3> 

CALL  SUB22CMT0T#M1#NB# DERI# DER3# RO#NESS>  

PR1NT3 60#  DER3 

3 60  FORMATC  ■ • THE  ANSWER  TO  FORMULA  £2  I S •#  DEC#  145 

PRI  NT  3 5 2#  NBS  Q 

CALL  VAR2C  DER3#  P#  VARN#  PO#  MTQT5  

3 52  FORMATC  »OTHE  EXPECTED  VALUE  OF  THE  DQUARE  0F_ THE  BACKORDER  LEVEL1 

1 St* # D20# 14?  

PR2NT41 0# VARN 

4 1 0 FORMATC  * OTHE  VARIANCE  0F_  THE  EXPECTED  NUMBER  OF  UNITS  IN  RESUPPLY 

1IS»»#D20«14> 

DER3-0.0 ___  

CALL  VAR2C  PER 3# P»  RE S PS 8# PO#NTOT> 

PRINT420#  RESPSQ  

420  FORMATC^GTHE  EXPECTED  VALUE  OF  THE  SQUARE  OF  THE  NUMBER  OF  LNITSI 

IN  RESUPPLY  ISI •#  D20.1 4) 

CALL  SNAVEC NT0T#M1#  P#  PO#  SM5 
PR1NT430# SN 

A 


4 30  FORMATC  *OTHE  EXPECTED  STOCK  LEVEL  IS*#D2W1 4) 

CALL  VAR3CNT0T»M1#  P*  PO# VARS#  SN5  

_____ PRI  NT  44  0#  VAR  S 

440  FORMATC  #OTHE  VARIANCE  OF  THE  STOCK  LEVEL  IS  *#D20.145 

_U«0.0 _ 

CALL _VAR3CNTOT#Ml*P*PO»SSQ#U5  

4 25  FORMATC  * OTHE  EXPECTED VALUE  OF  THE  SQUARE  OF  THE  STOCK  LE  VEL  IStj 

1D2  0.  1 4 > 

PRI  NT 42 5*  SSQ 

CALL  MNAVECNTOTjMUPjPOjMN) 

PR1NT435#MN 

4 35 FORMATC  • OTHE  EXPECTED  NUMBER  IN  SERVI CE  I S *,  D20«  1 45 

CALL  VAR4C NT0T>M1#MN» F>  VARM*  P05 

PRINT444#  VARM 

444  FORMATC  " OTHE  VARIANCE  OF  THE  EXPECTED  NUMBER  IN_  SERVICE  IS_*#D20i 

145 

CALL  VAR4C  NTQT#M  1 # U#  P#MNSQ#  P05  

PRINT446#MNSQ __  

4 46 FORMATC* OTHE  EXPECTED  VALUE  OF  THE  SQUARE  OF  THE  NUMBER  _IN  __SERViE 

1 JS«*#D20*14>  

STOP 

END 

SUBROUTINENINTENC  NTOT#M i#NB 1 * NB* PO 1#  PO# AN S 19* MO#  PA#  Rp5 

DOUBLE  PRECISION  NTOT#NBSQ#M  1#NB*NB  1#  PO# PO 1*  ANS  19#M0*  PAC  1 01  >»  RO 

ANS19-0.0 

CALL  VAR1C P01#Ml#NT0T* 1#NBSQ# PA# ANSI 95 

ANS19-U-C  l./CNT0T-MH-15  5*CCCM0*R05*P01/P05-NBSQ/NB5 


RETURN 


END 


SUBROUTINE  SUBS (XNTOT#  RO*  DER 1 * DER2 > 


C FINDS  THE  RESULTS  FOR  FORMULA  20* 


DOUBLE  PRECISION  NTOT* RO* DERI* DEH2C 101 > 


I TOT-NTOT 

DO  10  NRESUP— l® I TOT 


1 0 


DER2<  NRESUP >-NRESUP/RO- DER 1 /RO 


P0DER-DER1 /<  0*  O-RO  > 

RETURN 

END 

SUBROUTINE  SUB21C NTOT# P*  DER 1 ) 

C FINDS  THE  EXPECTED  NUMBER  OF  UNITSIN  RESUPPLY. 

DOUBLE  PRECISION  NTOT* PC  101 5* DERI  

DER1-G.0 

I TOT-NTOT 


7 


2350 


1 0 


DO  10  NRESUP- 1*1 TOT  

DERI -NRESUP* P< NRESUP >♦ DER 1 
RETURN 


END 


SUBROUTINE  SUB22<NT0T*M1*NB*  DERI*  DER2*  RO*NBSQ> 

C FINDS  THE  RESULTS  FOR  FORMULA  22* 


DOUBLE  PRECISION  NTOT*M 1 *NB* DER 1* DER2* RO*NBSQ 


SUBROUTINE  VAR3<NTOT#M 1# P# PO#  VAR#  SN> 


DOUBLE  PRECISION  NTOT#Mi#  P<  1 01  >#  PO#  VAR#  SN 


ITEST-NT0T-M1 


VAR-P0*< NTOT-M 1 - SN >**2 
DO  10  NRESUP— 1# I TEST 


1 0 


VAR- VAR*  P<  NRESUP>*<  NTOT-M 1 -NRESUP- SN >**  2 


ITOT-NTOT 


INI TAL-I TEST* 1 


30 


DO  30  NRESUP- INI TAL# I TOT 
VAR- VAR*  P<  NRSSUP5* SN**  2 

RETURN 

END 


SUBROUTl NE  SNAVEC  NTOT# M 1 # P#  PO#  SN5 

DOUBLE  PRECISI ON  NTQT#M1#P<_1 01)#P0#  SN 

ITEST-NTOT-Ml 

SN-PO*<  NTOT-M  15 

I F<  1TEST.LT.  15Q0T0  20 

DO  JL  0 NRESUP- 1 # I TEST 

SN— SN*P<NRESUP>*<  NTOT-M  1 -NRESUP/ 

1 0 CONTINUE 

2 0 RETURN 

END 

SUBROUTINE  MNAVEC NTOT#M  1 » P# PO# MN>_ 

DOUBLE  PRECISION  NTOT#M 1# P< 1 01 ># PO#MN 

INI TAL-NTOT-M 1* 1 . . 

MN-0 


I F<  INI TAL.LE. 0)MN-P0*NT0T 


ITOT-NTOT 

I F< INI TAL  »LE. 03 INI TAL-1 

DO  lONRESUP-lNI TAL* I TOT 

1 0 MN-MN*<NTOT-NRESUP)*P< NRESUP? 

MN-MN*M1»P0 

I TEST-  INI  TAL-1 

miTEST.LT.  1 ? GOTO  30 

DO  20  NRESUP- 1*1 TEST 

2 0 MN»MN*M  1 * P<  MRESUP  5 

3 0 RETURN  _____ 

END 

SUBROUTINE  VAR4<NT0T*M1*HN*P* VAR*  P05 

DOUBLE  PRECISION  NT0T*M1*MN*P<  1 01  3*  VAR*  PO 

1NITAL-NT0T-M1 

I TEST- INI TAL* 1 

VAR-PO  _ 

DO  20  NRESUP- 1*1  TEST 

2 0 VAR- VAR*  P<  NRESUP? 

VAR- VAR* < M 1 -HN  ?»»  2 

IF<INITAL.EQ.O>VAR-VAR*PO*<NTOT-HN?»*£ 

I F< INI TAL • EQ • 0 5 IN I TAL  » 1 

ITOT-NTOT  

DO  10  NRESUP- INI TAL* I TOT 

10  VAR-VAR*P<NRESUP>*CNT0T-NRESUP-MN)**2 

RETURN 


SUBROUTINE  FLOW<X#vJ»  IX# IY>  

C CHECKS  FOR  OVERFLOW  AND  UNDERFLOW  — PRINTS  WARNING  MESSAGE  1 F ONE  I SF 

0 UND* 

DOUBLE  PRECISION  X 

1F<X#LT»»174» 9> IX-IX+ 1 

1F<X«(3T.  174.9)IY»IY+  1 

IF<X*LT*»  17  4*9  > PRINTS#  J 

IFCX#LT*-174.9>x—  174*9  

5 FORMAT<  ■ »»»»»UNDERFLOW»»»»»  _ 

IFCX.GT* 174*9  > PRINT  1 0#  J 

10 FORMAT* • • *****OVERFLOW*****  J-«#I3)  

IF<X*QT*  1 74*9  ?X»  1 74*9 

RETURN  _ 

END 

SUBROUTINE  VAR1C PO#M 1#NT0T# VAR#  P# NB9  

DOUBLE  PRECISION  PC  1 01 /# PQ#M1#MT0T# VAR#NB  

ITOT-NTOT 

INITAL-NT0T-M1+1 

VAR-O.O  _ 

IF<INITAL.EQ.0>VAR-VAR+P0*CH1-NT0T-NB>**2 


I 


I F<  1NITAL»EQ.  0 9INITAL- 1 

DO  lONRESUP-lNlTAL#  I TOT  

1 0 VAR-VAR*PCNRESUP>*CM 1 -NT0T+NRESUP-ME>**2 

ITEST-NT0T-M1 


DO  20  NRESUP" 1# I TEST 


20 


VAR- VAR*NB**  2*  PC  NRE  SU  P ) 


VAR- VAR*  P0»NB»»  2 

RETURN _ 

END 

SUBROUTINE  VAR2C DER 1 * P#  VAR# PO#NTOT> 

DOUBLE  PRECISI ON  PC 10 l ># DERI#  VAR# PO#MTOT 

ITOT-NTOT 

VAR-  P0»  DER  1**2 

DO _1  0 NRESUP- 1 # I TOT ... 

1 0 VAR— VAR* PC  NRESU  P>fC  NRESUP-  DER  1 2 _ 

RETURN 

BID 

SUBROUTINE  PVALC  P#NTpT# HO#  RO#  P0> 

DOUBLE  PRECISION  ELOQ# DEXP# P0#M0#R0#X#C#B#A# PC  101  ?#NIQT 

ITEST-NTOT  

PO-DEXP<-MO>pRO)  

X-PO 

DO  10  NRESUP- 1, 1 TEST  

B-0«  0 

DO  30  J- 1#NRESUP  

A-J ' _ _ __  _ _ 

30  B-B*DLOQCA> 

— - --  ‘ 

PCNRESUP?-<MO»BO>»»NRESUP»  DEX  PC  -MO*RO  > / DEX  PC  B> 

X-X*PCNRESUP> 

PRINT20# NRESUP#  PCNRESUP> 
i 0 CONTINUE 


A-  25 


20 


40 


FORMATC*  • P<*»13»*>  - *#D20*14> 
PRINT40#X 


I 


FORMA  T<  • OTHE  TOTAL  OF  THE  POSSION  PROBABILITIES  IS*#D20«14_5 


RETURN 


END 


SUBROUTINE  NBSTXC NB# STOCK#NTOT*MO#Ml#RO»  P0> 


DOUBLE  PRECISION  EL0G»  DEXP»NTOT#Mlj»NB#  STOCX#  D»C*M0»  RO#  P0_ 
M-NT0T-M1+1  . 


I TEST- 100 


NB-STOCK-O. 0 
I FCM#LE#  0)NB— PO 


I F<M«LE»  0>M- 1 

DO  50  NRESUP-M# I TEST  _ 

C-0.0 

I F<NRESUP*LE*  1 >QOTO  65 


£ O 


60 

65 


DO  60  J-2*NRESUP 

D-J 

C-C+DLOGC  D) 


IF<NRESUP*ELOG<MO*RO>-C-MO*RO#LT.f  IGObOTOl  1 1 


5 0 NB-NB-KMl-NT0T+NRESUP)<»»DEXP<NRESUP»£L0G<M0»R0)-C3*DEXFt-M0*Rp) 

1 1 1 CONTINUE  


PRINT 100#NB 


100  FORMAT<  • OTHERE  ARE*# D20*  1 4»  •.  BACKORDERS* >_ 


RETURN 


END 

S ENTRY 

20  *2 


15 


17 


/LOGOFF 


CONSERVATIVE  SINGLE  SERVER 


INVENTORY  MODEL 

* (FINITE  POPULATION) 

PROGRAM  LISTING 


) 


1 


A-  2 7 


I 

. _ 

/LOGON  "NC"  

/CO  RUN.WATFIV 

S JOB  (L3NTJQQ5 __  KOVACS  

DOUBLE  PRECISION  PC  1 01  >#  PA<  1 G 1 >#  DER2C  1 01  >#NB#NB  1#  P0»  Pp  U2#NT0T#  Hi 

1M  1#M0#  AN  S3  6#  DERI#  DER3#  VAR#  VARA#  NB  SO#  VARN#  RESFSO#  VARS#  U#  SSQ#MNSQ#  J 

2ANS19 

Z-G#0 

PRINTS 

READ#  X 1 # X 2#  K 3#  X 4 

5 FORMATC • 1 •#  5X#  CONSERVATIVE  SINGLE  SERVER  INVENTORY  MODEL 

READ#  MIN#  MAX#RO#MO#M  1 

PRINT  1 0#MIN#  MAX#  RO#MO#M  1 

1 0 FORMATC  •CMIN-#  I4»  AX,  •MAX"*#  1 4#  4X#  «R0- •#  D2C.  14#  4X#  *M0"*#  D20#  1 4#  4> 

1 "Ml"*#  D20#14> 

DO  1 000  I TOT-MI  N#  MAX  

NTOT-I TOT  

IF<X4#EQ#  1>  PRINT}  50 

IX-IY-C  

I FCMO*RO«LE«  1 • OICALL  FVALC P#NTOT#MO#RO# F0#K2> 

IFCMO+RO.OT . 1.03CALL  PVAL1AC  P#NT0T#M0#R0#  P0#IX#IY#K2> 

C IX  AND  IY  ARE  TOTALS  FOR  UNDER  AND  OVERFLOWS  IN  THE  PROGRAM# 

1 50  FORMATC  •OPCN)  IS  THE  PROBABILITY  THAT  N ITEMS  ARE  IN  RESUPPLY  AST 

II ME  NEARS  I NF I NITY«») 

C T IS  THE  TOTAL  OF  THE  PROBABILITIES  — SHOULD  EQUAL  1 • OCC#  

T-0#  

DO  100  I- 1# I TOT 
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1 00  T-T+PC2  5 

T»T»PQ 

IF<K2#EQ.15PR2NT16C.T 

1 60 F0RMATC  *0THE  TOTAL  OF  ALL  THE  P VALUES  IS  '# D2C. 1 45  

CALL  SUB17(NT0T.M1.P.X1 5 

CALL  SUB13<NTQT.M1.M0.RQ.NB.X15 

I6-NT0T-M1+J 

2- 0.0 

DO  300  2 7"  I 6.  2 TOT 

3 00  Z-Z+PC  175 j 

C ZIS  THE  TOTAL  OF _ ALL  THE  PR0EAB2LI TIES  FROM  I-NTOT-M  1?_  1 „ TO  2 “N TOT. 

2FCX1  .EQ.  1 5PR1NT301.Z 

3 01 FORMATC/OTHE  PROBABILITY  FOR  STOCKOUT  2 S". 020.1 45 

CALL  STOCXCMO. RO.M  l.NTOT.X 1 5 

2X»IY«0 

KS-0  

IF<MO*RO.LE. U 5 CALL  PVALC  PA.NTOT*  l.MO.RO.  P01.K55 

I FXMO* RO.OT. 1.05 CALL  _PVAL1A<  PA. NT QT*  1 . M O.RO.  PO  1.IX.  2 Y.  K 5 5 

279-2TOT-M  

PRINT220.  IX.  2 Y 

2 20  FORMATC  * CTHERE  WERE  *.2  3.*  UNDERFLOWS  AND  *.2  3.  * OVERFLOWS*? 

t-o*o 

DO  102  !■  1.279 

1 02  T*T»PA<  2 5 

CALL  SUB1 SCNTOT*  1 .M 1 . MO.  RO.  NB 1 . K5  5 

CALL  SUB36(M0.R0. PO. P01.ANS365 

2 F<K4.EQ.  1 5PRINT400. 1 TOT.  2 TOT* 1. 2 T0T.ANS36  
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I 


L 


■ 1 

1 

4 00  FORMAT!  • CNBC  * . I 3.  * ) -NB!*.I3.*)  / NB!*.I3.*)  « *.D20.14) \ 

IFCX4.EQ.1)PRINT304.«-!NB-NB1>/NB  

3 04  FORMAT!  *OTHE  LEFT  SIDE  OF  THE  EQUATION  GIVES  *,  DSC.  14) 

CALL  SUB 33! NTOT.  P.  DER 1 ) 

2 FCK3.EQ.1J)  PRINT  31 0.  DER  1 

3 10  FORMAT!  *OTHE  EXPECTED  NUMBER  0FUN2  TS  IN  RESUPPLY  I S*.  DSC.  14? 

IFCK4. EQ.  DPRINT320.  - DER  1 /RO 

3 20 FORMAT!  * GT3E  RESULT  OF  FORMULA  33  IS  •.D20.14) ) 

CALL  SUB37CNT0T.R0. DERI. DER2) 

‘ 1 L 

2FCK4.EQ.  1 )PRINT333 

333 FORMAT!  •OTHE  ANSWERS  FOR  FORMULA  37*  •) | 

IFCX4.EQ.1 )PRINT330. 0. -DER 1 /RO 

IFCX4.EQ. 1 ? PR2NT330. ! J5. DERS! J5 )» J5- 1. 2 TOT) 

330  FORMAT!  * PC  *.  1 3.  * ) GIVES  *#020.14) j 

CALL  VAR 1 ! P0.M1. NTOT. VAR.  P.MB)  

I FCX3.EQ. 1 )PRINT350. 1 TOT.  VAR  ; 

3 50 FOHMAT! • THE  VARIENCE  OF  NB  FOR  NTOT-*. I 3.  • 1 S *.  D2G. 14) 

DER3-0.0 

CALL  VARl!  PO.M  l.NTOT.NBSQ.  P»  DER3)  

CALL  SUB39CNT0T.M l.NB. DERI#  DER3. RO.NESQ)  

IFCK4.EQ.DPRINT360.DER3 

3 60  FORMAT!*  OTHE  ANSWER  TO  FORMULA  39  IS  *.D20.  14) i 

2FCK3.EQ. 1 )PRINT352.NBSQ 


CALL  VAR2! DERI. P.VARN.PO. NTOT) 


E- 


4 10  FORMATC  *OTHE  VARIANCE  OF  THE  EXPECTED  NUMBER  OF  UNI TS  IN  RESUPFLT 

IIS*  •#D20*145 

DER3-C* 

CALL  VAR2<DER3#P#RESFSQ#P0»NT0T5  

I FC  X 3 • E Q • 1 5 PR  I NT 42 0*  RE  S P S Q 

420  FORMATC  *OTHE  EXPECTED  VALUE  OFTHE  SPUARE  OF  THE  NUMBER  OF  UNI  TSI 

IN  RESUPPLY  ISf»D20»145  

CALL  SNAVECNT0T»HUP»PQ»SN5  

I F<X3»EQ»  1 5FRINT430*  SN 

430  FORMATC  * OTHE  EXPECTED  STOCK  LEVEL  IS  •#  DSC*  145  

CALL  VAR3<NT0T#M1,  P#P0#VARS,  SN5  

I F<  K3>  EQ»  1 ) PR  I NT 44 0«VAR  S 

440  FORMAT<  • OTHE  VARIANCE  OF  THE  STOCK  LEVEL  IS  •«  020*145 

u-q. 

CALL  VAR3C NTOT#M 1# P# PO# SSQ#U> 

425  FORMATC  * OTHE  EXPECTED  VALUE  OF  THE  SQUARE  OF  THE  STOCK  LEVEL_I  SI# 

1D20.145 

I F<K3«EQ*  1 5PRINT425*  SSQ 

CALL  MNAVEC NTOT# M 1#  P# PO#MN  5 

IFCK3#EQ.15PRINT435#MN  

4 35  FORMATC  ■OTHE  EXPECTED  NUMBER  IN  SERVICE  IS  *jD20>145 

CALL  V AR4C  NTO  T«  M 1 # MN#  P#  V ARM# P05 

IFCK3.EQ. 1 5 PRINT 44 4# VARM 

44.4 FORMATC  * OTHE  VARIANCE  OF  THE  EXPECTED  NUMBER  IN  SERVICE  I Si  *#  PEG* 

145 

CALL  VAR4<NTQT*M1»U»P»MNSQ#P05 
IFCiC3#EQ.15PRlNT446#MNSQ 
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446  FORMATC  • OTHE  EXPECTED  VALUE  OF  THE  SQUARE  OF  THE  NUMBER  IN  SERVIB 


I I SI  *#D20#I43 

1 OOO CONTINUE 

sjrop 

END 

SUBROUTINE  FVALC  P#NTOT#MO#  RO#PO#X23 

DOUBLE  PRECISION  P<  10 1 3* N TOT# MO# RO#  PQ#Q#A# B#  C#  D# E 

DOUBLE  PREC I SI ON  CL 00#  DEX P 

ITEST«NT0T*1 

no  5 JK" 1#  2 TEST 

NRESUP-JK- 1 

CALL  QMQR0< MO#  RO#  03 

M»M0 

CALL  FACTCM# C 3 

A—M0»CLQq<M03+C 

B-NTOT*  EL  OQC  M 0*_R0  3 

A“A+B+ti  

B-0>0 

I-NTOT-MO+1  

DO  10  M8-1#I 

J-M8-1 

10  B-B*DEXP<J*CL0tKM0*R03  3 
A— C CLOCK  DEXFC  A3*B33 


ix-ljr-q 

1 FCNRESUP* EQ#  0*  03 CALL  FuOVK  A#NRESUP#  IX#  IY3 
IF<NRESUP#EQ*0*03F0-DEXP<  A3 


I F(MRESUP*EQ# 0*  05 GOTO  5 

2 F<  NRESUP.LE.NT0T-M03P<  NRESUP3-DEXP<  ( NRESUF*  ELOM  MO*  RO  > )♦  A) 
D-<NT0T-M03*EL0G<M03 


MSMNTOT-NRESUP  

G-0.0 

CALL  FACT<M9»6> 

E*MR£SUP»CLOeC R03 

IFCNRESUP.  GT.NT0T-M03  P<  NRESUF3-DEXP<  C-G+D+E+A  3 

2FCK  2#EQ  *13  PR  I NT101#  MRESUPj  PC  N RESUP  3 

5 CONTINUE 

101  FORMATC*  P<»*23*»3  - •#D20*1A3 

RETURN 

HUD _ _ 

SUBROUTINE  SUB 1 7< N TO T# M 1 * F#X13  

DOUBLE  PRECISION  NTOT*MUF<  1013»NS*X  

NB-0*0_ _ ^ 

M-Ml  __  ___  _ 

DO  10  J-UM 

X-J-l 

I-NTOT-X  

NB«NB*<M1-X3»P< 13  

10 CONTINUE 

I F<  K 1 • EQ*  1 3 PRI  NT  AG#  NB 

AO  FORMATC  * OFORMULA  17  01 VES»* D2C. 1 A* * BACK  ORDERS  ^3 
RETURN 
END 
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SUBROUTI NE  SUB3 6< MO*  RO*  PO*  PO l * AN  S3  6> 


DOUBLE  PRECISION  MO*RO* PO* P01*ANS36 


ANS36-1*0-M0*R0*P01/P0 


RETURN 


SUBROUTINE  SUB37CNT0T* RO*  DERI * DER2 > 

DOUBLE  PRECISION  NTOT* RO* DERI* DER2C 101) 


ITOT-NTOT 


DO  10  NRESUP-  1*2  TOT 


DER2<NRESUP)-NRESUP/R0-DER1/R0 


PODER-DER 1 /<  0.  0- RO 5 


RETURN 


SUBROUTINE  SUE33<MT0T* P» DERI > 

DOUBLE  PRECISION  P< 101 >*NTOT* DERI 


DER1-0.C 


ITOT-NTOT 


DO  10  NRESUP- 1* I TOT 

DER 1 -NRESUP1*  P<  NRESUP  J-t-  DER 1 


RETURN 


SUBROUTINE  SUB39C NTOT* M 1*NB* DERI* DER2* RO*NBSQ> 


DOUBLE  PRECISION  NTOT*M 1*NB*  DER 1 * DER2* RO* NB SQ 
DER2-NB/R0*<<NT0T-M1)-DER1*NBSQ/NB> 


RETURN 
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SUBROUTINE  VAR3<NT0T#Ml#P#P0#  VAR#  SN> 


k 


DOUBLE  PRECISION  NTOT#M  1 # PC  lOl  3# PO# VAR#  SN 

ITEST-NTOT-MJ 

VAR«P0*<NT0T-M1-SN3**2 _ 

DO  10  NRESUP* 1# I TEST 

VAR-VAR*P<NRESUP3*<M-M1-MRESUP-SN3**2 

ITOT-NTOT 

INI TAL-I TEST* 1 

DO  30  NRESUP-INI TAL# I TOT 

VAR-VAR-fPC  NRE  SUP3*  SNmm  2 

RETURN 

BID 

SUBROUTINE  SNAVEC NTOT#M 1# P# PO# SN3 

DOUBLE  PRECISION  NT0T#M1#P<  1QO#PO#SN 

ITEST-NTOT-Ml  

SN-PO*C  NTOT-M 1 3 

ifcitest.lt. n e o to  20 

DO  lO  NRESUP-1# ITEST 

SN-SN4PC  NRESUP3»< NTOT-H 1 -NRESUP3 

1 0 CONTINUE 

2 0 RETURN 

END 

SUBROUTINE  MNAVE<NT0T#M1#P#P0#MN3 

DOUBLE  PREC1 SI ON  NT0T#M l#  PC  1 0 1 3# PO# MN 

INI  TAL -NTOT-M  14  I 

MN-0 

I FC  I N I TAL  • E Q.  0 3MN-  PO*  C N TO  T 3 


ITOT-NTOT 

IFCIN  ITAL.EC.C?  IN  I TAL,  - 1 
DO  1 0 NRESUP-INI  TAL* I TOT 


10 MN-MN+C  NTOT-NRESUP?>#PC  NRESUP? 

MN-MN+M  1*P0 

_ I TEST* INI  TAL- 1 

IFCITEST.LT* 1 )QOTO  30  

DO  20  NRESUP* 1#I  "EST 

SO MN-MN+Mljn  PC  NRESUP?  

3 0 RETURN 

END 

SUBROUTINE  VAR4CNT0T.M  l.MN#  P* VAR*  PO? 

DOUBLE  PRECISION  NTQT#M1#MN# PC  1G1?#VAR#P0 

INITAL-NT0T-M1+1 

I TEST-INI  TAL- 1 

VAR-PO 

DO  20  NRESUP* 1*1 TEST 

20  VAR-VAR^  PC  NRESUP? 

VAR-VAR*CM1-MN?**2  _ 

IFC INITAL.EQ.O?VAR*VAR»PO»C NTOT-MN  ?**2 

I FC  INI  TAL.  EQ.  C?  INI  TAL*  1 

ITOT-NTOT  

DO  10  NRESUP* IN I TAL# I TOT 

1 0 VAR* VAR#  PC NRESU  P?»  C N TOT-N  RE  SUP-MN?** 2 

RETURN 


A-  36 





END 


SUBROUTINE  FLOV<X,J# IX# I Y> 

DOUBLE  PRECISION  X 

IF<X.LT«-174«  9 5IX- 1 X+l 


I F<  X*  eT.174.9  >X-1 7 4.  9 

RETURN 

END 

SUBROUTINE  VARUPQjHUNTOT»UAR»P»NB? 

DOUBLE  PRECISION  P5  101  >,PO»MUNTOT#  VAR*NB 

ITOT-NTQT 

INITAL-NTOT-Ml+1 

VAR- 0*0 

I F<  INITAL.EQ.  0 )V AR-VA R+FO-  <M  1 -NT0T-NB5— £ 

IFClxVI  TAL  «EQ«  OINI TAL- 1 

DO  10  NRESUP— IN1TAL*  I TOT 

1 0 VAR»VAR*P<NRESUP>*CMi-NTOT*NRESUP-NB>**2 

ITEST-NTOT-Ml 

DO  20  NRESUP-UITEST 

8 0 VAR-VAR+N B**2* P< NRESU P > 


VAR-VAR* P0«NB««  2 

RETURN 

END 

fc 

A-  3 7 

. 

L 

2 

SUBROUTINE  VAR2C DER 1 # P# VAR# PO#NTOT> 


DOUBLE  PREC2 SION  P<  101  >»  DERI# VAR# PO#NTOT 

ITOT*NTQT 

VAR- PO* DERI** 2 

DO  10  N RE  SUP- 1 # 2 TOT 

1 0 VAR-VAR+P<NRESUP)*CNRESU P- DERI)** 2 

RETURN 

END 

SUBROUTINE  SUB19<NB#MO#RO#NTOT#M1#ANS19)  

DOUBLE  PRECISION  M0#R0#NB#M1#NT0T#ANS19#B#B» A 

DOUBLE  FRECI SION  0.0(8# DEXP 

CALL  QM0B0<  M0» 1.0#  RO#  fa ) 

G-DEXF<G> 

M-MO 

IFCMO-1.0.LE.O.O)Q*Q.  0 

CALL  F ACT  CM#  A > 

B-MN-NB-1  • C/RO 

IFCB.LT.  1 E-12.  AND. B • QT.0.03  B-  0.0 

A*<  (1  • 0-M0*R0)*DEXP<  -MO*£LOb(MO)*A)*Q*E)*<M  l-CMO*  RO)*(M  1-NB)  )«DES 
U CU  G-M03*EL0G<M0*R0)  )♦<<  1 . O-DEXPfC  M1-M0)*EL0C«<M0*R0>  ) >*DEXK  -M  1« 

2L08<  M0*R0))/C  1 . 0-M0*R0> >-<M 1-MO  )*DEXP< -MO*ELOB<MO*RO  > ) 

A MS  19*  C DL  OOC  A)  »£L  OMfJB)  )/<-EL  OCKMO*  R 0 » 

RETURN 

END 

SUBROUTINE  FACTC M#A) 

DOUBLE  PRECISION  A# B 


DOUBLE  PRECISION  CLOU 


A-O.C 


k 


I FCM.LE.  1 >OOTO  20  

DO  10  J-2«M 

B-J 

1 0 A-A+OLO0CB> 

2<^ RETURN 

HID 

SUBROUTINE  0MOROCMO.  RO.  A>  

DOUBLE  PRECISION  MO»RO»A»  B.  EL  00.  DEXP 

INTEGER  X 

M-MO 

A-O.C 


A-ELO0<AJ 


CALL  FACTCM»A) 


B-Ml  - NB-  U 0/R  0 

IF<B*LT*1E-12*AND*B*QT*  0*  0>8**C«  G 

A-<  CM3-<M0*R0>*v-.0-NB>>*DEXP<<  1 • 0-M03*  DL0Q<M0*R0>  ) >♦(  DEXP<-M0*DLQ 

1<M0J*A>*G*C  W0-M0*R03*B> 

ANS-C  EL0G< NB)-DL06CA>>/<  EL  0Q<  M0*R0>3 

PRINT 10* AMS 

10  FORMATC  * OFORMULA  21  8IVES  *»D20«14> 

2 0 RETURN 

END t 

SUBROUTINE  PVAL1 ACP# NTOT#HO»BO»PQ» IX#IY#K23 

DOUBLE  PRECI SION  PC  1 01  >#NTOT*MO«RO» PO# &#A*B#C# D«E*F 

DOUBLE  PRECISION  BLOtJjDEXP  

I TEST ■ N TO T+  l 

DO  5 JK« 1*1 TEST  

NRESUP-JK- 1 


CALL  GMORO<MO#  R0»  t»> 


M-MO 

CALL  FACTCM#C> 

A«-a0»  EL  OQC  MO  >♦ C» ft 

B-0*0  

I-NTOT-MO+l 

IF< I.LE.O)QOTO  20 
DO  10  M3«l»2 
J-M8- 1 


B«B* DEX  P< ( J-NTOT >* ELOGC MO*  RO  > > 


A— <EL0O<DEXP<  A5+B55  

1 FC  NRESUP#  EQ#  0*  0>CALL  FLOWC  A#  NRESUP#  IX#  I Y } 

IFCNRESU  P .E  Q •0)P0-  DEXPtC  -N  TO  T*  EL  0 G<  M0*R0  5*A  5 

2 FC  N R E SUP  • E Q • 0 • AND*  K 2 • E Q • 1 5 PRINT1C1#NRESUP#  PO 

I FC  NRESUP#  E Q • 0#  0 5 Q0T0  5 

IFCNRESUP«LE#NT0T-M0?PcNBESUP5«DEXPCCNRESUP-NT0T?»EL0toCM0»B0?+A> 

D— M0»EL0GCM05 

M9-NT0T-NRESUP 

F-0.0 

CALL  FACTCM9#  F> 

E-C  NRESUP-NT0T5*  ELOGC  RO) 

I FC  NRESUP#GT#NTOT-MO>  PC  NRESU  P>- DEXPC  C- F*  D*E* A} 

IFCK2#EQ#1>PRINT  1 0 UNRESUP#  PC  NRESUP? 

5 CONTINUE 

1 01 FORMATC  * PC  *# 1 3#  * > - *#D20# 145 

RETURN 

END 

SUBROUTINE  SUB13CNT0T# Ml # MO#  R 0# NE#K15 

DOUELE  PRECISION  NTOT#MG#M  1#  RO#NB*  Q#  GAMMA#  BETA#  A#  E#  C#  D#E»  F#G#Gl#H 

1 R#  S#T#U#  DLOG#  DEXP 

M-MO 

e«o#o 

IFCMO.LE# 15Q0T0  5 

CALL  QM0~9C  M0»  1 » RQ#  8 5 __  

G-DEXPCU) 


A-4 1 





vp-V 


r 


5 CALL  FACT<M*U> 

A-<M1-1*  0 /RO  ><•  DEX  P<-M  O-EL  OlKM  0 )♦  U >*  G 

B-DEXP<  CLOCK M1)>K  1 • G-MO>*ELOG<MO*RO>  > 

GAMMA-MI -MO 

BETA- 1 > 0-M0» RO 

C-0.0 

M9-GAMMA- 1 

IF<6AMMA*LT*2*G>Q0T0  20  

DO  10  J-  UM9 

K-J-l 

CALL  FACTCM9* 2 

CALL  FACT< K» E > 

J5-6AMMA-K-1  

CALL  FACTC J5#  F? 

H-K+2 

1 0 C-C+C- 1 3— K*6ETA**K*  DEX  PC  D-E-  F-OLOG<H3  ) 

20  R-0*0 

I Ft  BETA*  EQ*  0*0*  OR* GAMMA*  E Q • 0 • 03 GOTO  25 

R-C - 1 )-*<M9*£3-EETA-*M9- GAMMA  

25  2F<M1*EQ*1*C>Q-C«R 

1FCM1  »ME*  1 *0JQ-DEXP<<  1*  0-M1  >*EL0G<M0»R0>  C+R3 

J26-NT0T-M0* 1 

S-0*0 

DO  30  J27-UJ26 

1 

30 S-S+C  DEXP<<  J-MTOT?*  CLOCK M0*R0>  > > 

CALL  GM0R0CM0* RO*  G 1 3 
T-DEXP<-MO*ELOGCMO>+U+Gl) 
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■ ■ 


1 


r 1 

NB-DEXP<EL0B<A4>B*Q>-EL0e<S«>T3>  

I FCX 1 • EQ#  1 5PRINT50#NB  

IFCKUEQ.  1 5PRINT70#  + DEXP<  ELOC><  A*B*Q>-ELOtK  T)> j 

FORMATC « 0 FORMULA  13  GIVES »» D20* 1 • BACKORDER S * ? 

FORMA T<  * OFOR  N- INFINITY*  NE  -•#D20#14>  

RETURN 

END 

— 

SUBROUTINE  STOCKCMO#RO*Ml#NTOT»Xl ? 

DOUBLE  PRECISION  HO» RO*M l#NTOT# Q# A, B, C# D# E#  F#  STOX 

DOUBL E PRECISION  J3LOO* DEX P 

M— MO 

CALL  FACT<M»A> 

CALL  BMOROC  MO#  RQ*  O > 

_ B-DEXP<-MO»ELOQ<HO?-»A+Ci> 

C-0.0 

I -Ml -MO 

IF< I • L T • 1 5 6 0 TO  20 

DO  10  Nf-  1*  I 

J-N-l 

I FCM  1-J.  EQ.UOC-C*  1 

I F<M 1 - J, EQ.l .05Q0T0  10 

C-C+DEXP<  < J-M  !♦  1 >»EL0OCM0»B0>  ) 

10 CONTINUE 

2 0 D—  0*  0 

I J-NTOT-MO* 1 
DO  30  N-1#IJ  __ 

J-N-l 


7 0 
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IF<MO*RO»  ST.  l • 0>ST0K-DEXP<  £L06<  E*«  B+O5-EL0QC  D+B*E5> 

IFCKUEO.l  > PRI NT5  0»  STOX 

5 0 FORMAT<  •01HE  STOCKOUT  PROBABILITY  IS«»D20*14> 


RETURN 


PROBLEM  1 
PROGRAM  LISTING 


A 
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OLURLE  PRECISION  N,NBO,U*L 

0 i U B L E . .PREC  I 5 I JNi  A>  6 ,C,  D,^;il  >ku2.>  RD  3 > R r;  4,  lo,  Rg,  A2 1 , f ■ r.  r T > 
pN_l_>  £K  2 *fH  3 >_  P (1 .0 1 ) > P □ > 3 / Q It  0 2 ^ «3  /q  4 ,QJ  ,06,  3 7,0  8 
INTEGER  DNUDN2 
READ  R0,M1,M0,NB0,NMAX,A21 

NM 1 11=  1 

N=NMAX/2 
IF(N.LT.l)  J*I 

A,a,C,0,R01,RC 2,  RG3jRQ4  = 0.G 

5 IFIN.EQ.  UGQTQ  13 

CALL  FIMD(N-1/RQ1^MI>Mo,HBO/A21) 

A- ( N-l ) »( 1 +R0/RQ1 ) 

10  CA1.L  FINOCN^R02jvH>MO^WBO#A21) 

3=N* ( 1+R0/RQ2  ) 

I_F  (N. EO. UMAX)  GO  TO  _30 

CALL  FlN0(.'J+l,Ra3,Ml,Ma,NBC,  A~2i  > 

C*<N+l)*< 1+R0/R03) 

30  _ IF(NMAX-1,J.EQ.N)GQTd  <*q 

CAUL  FINDM  + 2,  Ra4,Ml,Mn,NB0, A2i ) 

0=  (N+2)*<  l+RO/RO'4  ) 

40  FMjsC-B 

FNZ=8-A 

FN3=D-C 

g=N 

Qi  = B 
QZ=R02 

Qlztlzi 

Q4=A 

Q7=R01 

Q6=N+1 

Q7=C 

QS=K03 

IF(FNl.GE,0.0.AMD,FM2,LE.0.0)GdTC  100 
I F ( N. EQ .NMAX  . AND . FH2 . LE  .0 .0  )GU  ft)  100 
i F (N.EQ. 1.0. AND. FWl.GE. 0.0) GOTo  100 

_ Q»H+1 

Q1=C 

Q2=R03 

Q3 =N 

Q4»B 

Q5=RD2 

Q6»N+2 

Q7«D 

QR=R04 

IF (FN3.GC .0.0. AND.FN1 .LE.O.O)GuTO  100 
I F { FN1  . LE  .0 . 0 . AND ,NMAX-1~  0 . EQ.»'* ) GOTO  ICC 
Q*N+2 

Qjj*D 

Q2=RD4 
Q 3 = N -f  1 

Q*»C  ___  __  

Q5  »R03 

G6,Q7, 08=0,0 
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T 


I F ( HM  AX -2 .0  , E Q . M . AND.  FN3.LE,  0 . J ) G 0 T Q J , n u 

Q1*A 

Q2=RQ1 

Q3*Q4>05=0.0  ' ~ " 

<36  = 2 

Q7»B 

Q8=R02 

IF(N.EQ.2,0,AND.  Fri2.GE.0.0)CCTu  lQO 

1F(FN1.EG.Q)PRINT,  1 ER R flR  ft 1' 

IF (FN2.EQ. 0.0) PRINT* ' ERROR  #2' 

DNl  = l/a-(FN3/FNi)  ) 

DrJ2«I/.(l-JLFfU/FN2I  ) 

• I F ( DNI . LT  . ON 2 , AND . -ON  I . GT  . DM2  ) ..*N+ON  I 
lF(DiU.LT.DN2.AND.-Dhl.LE.Dti2)-«*N+l)M2 

IF(DN  2 s_L_E  . DN I.  AND  ^ DM  2 , GT . DM  1 )i<=N  + 0 N 2 

I F ( QN2  . UE  . DNI . AND  , -DM2  . Le.DNl)  N*N+DN1 
IF{N.LE»0)I=1 

I F(  R.  GIj  N M A XlRr  NM  AX 

PRINT*N*B*RQ2 
GOTO  5 

1 0 0 £EJL!Hl  Q0-»  01*  04*  05 

PR INT200*Q*Q1*  32 
PRINT200*Q6*Q7*Q8 

L*  ( 09 -Q2 ) V 02 ) / .(  (_a*i-Q8-0 21/02  ) - I I 

L«(  l-<Q2-Q5)/Q2)/( (Q*(02-C5)/Qc)-1) 

P°  INT21C-*l  *U 

PF  INTZ20*l .»&2,  J»Q2 

21 0 FORMAT { IX* 020. 14* » < ( R0/R01 ) < • * u20 . 1 A > 

220  FCKHAK  1X*D20.U*  ' <R0<  ' * D2C  , 14  J 

20  v FORMA  T ( 1 QN»  1 * D2D  .14*5  .X*  »Z='*D2u  . 14  *_/  * » RP  L = 1 *02  QjlU  ) 

PF  I NT 900 

900  FORMAT ('I') 

_ PAUSE  

STOP 

END 

SUBROUTINE  SUB 21  (?|TQT*P*DFR1) 

c finds  the  EXPECTED  NUMBER  of  UNITS  in  RESUPPLY. 

O.P.  NTOT*P(  Id)  *DER1 

•  DER1S0.0 

- ITQT*NTDf 

• DQ  10  NRESUP=1* I TGT 

'10  DER1»NRESUP»P (NRESUP )+QEPl 

' RETURN 

' END 

JUBRQUTIN  E HI  N E <_  'NT  0_T*M  0* P* PQ* RP* NB* Ml ) 

»C  FINDS  THE  PROBABILITIES  FOr  j*0  TO  j-fiTOT-MD  ( USES  NATURAL 
logsto  allow  for  LARGER  NTOT) . 

> DOUBLE  PRECISION  3 ( 1 C 1 ) * I ITCT*  Pu*  M(J,  RO*  A*  C * B*  D*  E*  F , G*  H*  7,  V*  NB  * M 

> ITEST  = NT0T-M0+1  

I11=ITEST 
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1 1 1 1 - .>.  JJ.UI.UJJUUmjUIUl|W»Jiii,  ,.v 


'iFdll.LT.  1)111  = 1 

DO  _10_  J1  = 1HTEST 

' ' JsJl-l" 

' Y.»  Z =0 . 0 

J f LIT E ST. IT , DSD T 0 2_i 

DC  20  M*1#ITeST 
NHESJP=Ml-l 

_ A_=  (NRESUP-J  ) *LGG(MP»RQ) 

C > D = 0 . G 

I F ( J . LE . D GQTQ  35 

DO  30  K = 2/  J 

B = K 

30  C=C+LOG(b) 

35  IF(NRESUP.UE.1)G3T0  <*5 

DC  40‘K»Z*NRESUP 
B = K 

<*0  D=D+LQG  ( B ) 

<♦5  1 E«A»C-D 

CALL  FLOW(E* J* IX,IY) 

IF  (E.LT.-174.99b)GDTQ 

Z*Z+EXP(E) 

■20  CONTINUE 

21  IFC  ITEST.LE.I)  IT,EST»1 

ITCJT*NT0T+1 

ITE3T=ITEST+1 

00  50  _N1  = I T E_S T j HOI 

NRESUP=N1-1 
00  = 0.0 

M = M 0 

I F ( M . L E . 1 > GQTC  65 
DO  60  K=2>i! 

B = K 

60  C*C+LOG<R> 

o5  I F ( J . LE . 1 ) GQTQ  75 

D0_70_  Kf2jJ 

B = K 

70  D=D+LDC(B) 

75 A ,y  F = 0 . 0 

IFJNRESUP.LE.l )GDTO  85 
DC  80  K-2j NRES JP 

B^K 

30  A = A+LQG ( B ) 

65  N*NTOT-NRESUP 

IF  ( H . L E jt  1 ) G C TO  95 

DO  90  K = 2 / N 
B*K 

90 _ F = F>LOG(B  ) 

95  G*(NTOf-MO-J)*LQG(NO) 

H*(NRESUP-J)*L3GCR0) 
E*C+D+C+M-A-F 
CALL  FLOUi  < £,  J,IX,IYI 
I F ( E . LT .-17A.996JG0T0  50 
Y= Y+EXP ( E ) 
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so  continue 

IF( J.EQ.O)PQ=L.O/ (Y>z  ) 

I F { J.EQ.0)GCTD  10 
P(J  )»i Y*Z) 

10  CONTINUE 

C A LL  T gN  (_N TpT t M 0 , p, R 0 ) 

CALL  S I XT E~'  1 ( NTOT , M l , P ; ;r p c ) 

RETURN 

END 

SUBROUTINE  TEN  ( NTOT/MC/ P* ro  ) 

C FINDS  THE  PROBABILITIES  FCR  j=NTDT-MG+1  Tq  J=imTL)T  (USES  NATURALLCGS 

TP  allow  for  larger  NTnri. 

D.P.P( 10l)#NT0T,Pa#H0*RD^A/B/CjD/E^ H#Zj Y 

IMTALsNTqT-MQ  + 1 

IF(  INITAL.LT.  D I.MITAL«1 

ITQT=NTGT 

DO  10  J»INXTAL>ITDT 
Y.»  Z = 0 ♦ 0 

DU  20  N 1 = 1 / J_N I T A L 

NRESU'P-Nl-1 

A^C/H,G=0.0 

I F ( J . L E » 1 ) GOTO  135 

DU  100  K*2*J 

S = K 

100  G=G+LOG<B) 

105  ' NsNTOT-j  ' 

IF(N.LE.1)GGTD  115 

DC  110  K=2j N 

B-K 

110  H = H + LOG l B ) 

115 If  itfRE S UP  , LE . 1 ? G3T0  35 

DC’  30  K = 2/f|RESJP 
B=K 

30  A = A -t- LU G ( E>  ) 

35  N = M U 

I F ( N . LE • 1 ) GCTO  45 

DL_  40  K = 2^,J 

B = K 

40  C=C-*-LOG(B) 


•O  Us  l liKtiUr-g  JVLLlulKriJ 

E = (NRESUP-(NTqT-'IO)  )*LUC(MQ) 

f«d+e+g+h-a-c 
call  flowcf, j, ix# r y j 

IF  (F.LT.-174.996)GDT0  20 

Y=Y+EXP(F) 

20  CUNTINUE  . . 

or  50  MRfcSUP=INITAL> ITOT 

A,OD/E/F=3.0 

IF( J .LF . 1 )G0TC65 

DU  60  K=2,J 
B=K 

oO  A = A + LQG ( B ) 


90  F=F+LQG(&) 

95 Ps<NRESUP-J)»L3G(RQ) 

B = D+A+C-~E-F 

CAwL  FLOWIB, J, IX, IY) 

IF<B.LT.-174.996)G0T0  50 

Z=Z+EXP(B) 

50  CONTINUE 

P ( J ) = I.0/( Y + Z ) 

10  CONTINUE 

RETURN 

END 

SUBROUT  I NE  S I X T E M ( NTQT, M 1 , P , NB> PC ) 

C Fir^OS  THE  EXPECTED  NUMBER  OF  BACKCkOERS. 

DOUBLE,  PR_EC  IS  ION  NTQT, Ml,  P(  101 ) jJJ&jPDj  3 

fTUT =NTOT 

■ t.B  = n. 

II>'ITAL=NT0T-M1  <-l 

■ * lFTTTnm.LE.O)NB  = PO 

* I F ( INITAL.LE.O) INITAL  = 1 

• 00_10  NR ESUP  = [NI TAL# I TOT  

B=(M1-(MT0T-NRESJP) )*P(NREbUP) 

10  NB=NB+& 

RETURN 

END  ‘ ~ 

SUBROUT  I N£  FLOW ( X, J, IX, I Y) 

C CHECKS FOR  OVERFLOW.  Af|D  UNDERFLOW  —PRINTS  WARNING  MESSAGE  IF 

JONE  IS  FOUND. 

D.P  .X 

I F (X.LT.-174.996 ) IX= IX+1 

I F(X.GT. 174.995) IYsIY+1 

IF (X.LT.- 174. $96) PR  I NTs,  J 
5 FORMAT  ( ' *****  UNDERFLOW*****  j=',I3) 

I F ( X . G T . 1,74.99.5 i PR  INTI 0,  j 

10  FORMAT  ( 1 *****0VERFLOW*****  J='M3) 

I F(X.GT. 174.995 )X= 174. 99 

C SETS  X EPUAL  T]  MAXIMUM  VALUE  WITH  NO  OVERFLOW  IF  AN  QVEPFLOW 
IS  FOUND  UN(  175)>10»*76)  . 

return 

END 

_ SUBROUTINE  NBSTK  < NB,  NTJ..T,  MC,  M 1,  RO  ) 
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D . P . NTCJT,  H 1,  NB  , S TOCK/  Dj  L j MC/ RU*  ?U 

m«ntot-mi+i 

_ ITEST=IOO 

NB»C.O 

IF<M.LE.u)NB  = E/<P(-MO*RU> 

I_F<_H.J,Ft01rl=l 

DC  50  HRESUP»M,ITEST 
C=0.0 

I F < NR E SUP  . LE « 1 ) GOTO  65 

DO  60  J = 2/i JRESUP 
D = J 

60  C=C>IQC<D) 

6 5 I F ( NRESUP*LDG  ( MO*RO  ) -C -Mr*RG . L T . - 100  ) GOTO  1 1 1 

50  NB  = NBi-(H1-NT3T*NRESUP  ) *EXP ( NR E jUP*10G ( MO*RU ) -C  )*CXP  (-nn*PQ  ) 

111  CONTINUE 

RETURN 

END 

SUBROUTINE  FIHD(NT0T>R;U,M1>MCjNB2A/A)  

DOUBLE  PRECISION  ROT^Bri/HTC'TTCuUNf  , MO,  PB  ( 101 ) , PQB*fiB2j 
PC  ( 10l')/P0C/>tB3/MiVDER5jDERlC/0R  \»  EPS  I LN*  A/ NB2A 
R0>R01=1 .0 

1032  IF(A.NF .O.OJCALL  N I NE ( NTDT- Ij Mu; P0, POB, RC 1 , NB2, Hi ) 

IF i June. o.o) call  nine<ntoTjMO*pop dc * r d 1 > n a 3 , m 1 ) 

IFIA.NE.O.OJCALL  SUB21(NT0T,PC*DER1C ) 

IF  c A .19. 0 .0  )_C_ A L L N 8 S TK(HB  3 * NT  QJ  Q ,_M1 , RQ1) 

IF(A.EQ.O.O)CALL  NBSTK  (NB2/MT[ji-I/MG/MDRCD 
IF(A.EQ.0.0)DER1C=M0*R01 

COUNT  »COUNT+JL 

DER5=(-DER1C/R31 )*<N83-Nb2) 

DRO= (NB2A-NB3 ) /DER5 

_RU_1  = R0 1 + DRQ_ : _ 

EPS  I LN=DRQ/ROl 

IF  ( EP  S ILN.LE.2E- 1A.AND.EPS  I LN.oE.-2E-D)  GOTO  1031 

I F ( COUNT .GE. 100) GOTO  1031 
'GOTO  1C32 
1021  RETURN 
END 
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/LOGON  "NC”  

/ DO  RUN  • V ATFI V — 

S JOB  <L8NTJQQ> KOVACS 

DOUBLE  PRECi  SI  ON  N#_H0P1 # R0P2#  ROA#  ROB#  ROC#  DROl  P#  DR02P#  A I#  A3 

DOUBLE  PRECI  SI  ON  A#  B#C#  D#  RO 1 # R02#  R03#  R04#Z  0#  RO#  A2 1#  M 1#  MO#  FN  1#jFN2# 

IFN3#P<  101  >#  PO#  0#  Q 1 #22#  03#  04#  05#  06#  07#  06 

INTEGER  DN1#EN2  

READC  5#»>  R0#M0#Z0 

DO  1000  M31«l#2 — 

M1"MQ - 

_ A21 -M31-1.0 

PRINT#R0#M0#Z0#A21 

NMIN-1 

NMAX-ZO 

IF<NMAX«EQ#ZO>NWAX*ZO-UO 

I -Z  0/2 

I FC  A2 1 • EQ#  0#  0>N-I 

IFCN#LT#1>M«I 

A-B-C-D-RO 1 -R02-R03-R04"  0.  0 

5 IF<N>EQ# 1 >QOTO  10 

R01«R0»<N-1?/C20-»<M-1  » 

IFCA21#EQ#0#0>CALL  NBSTKC  A#N*  1#M0#M  1#  RO  l > 

1 F< A2 1 • ME* 0 • 0 > CALL  NINECN-  1#M0# P# P0#R01# A#M  1 > 

10 R02«R0*N/<Z0»N  ) 

IF<A2WEQ.0#0>CALL  NBSTKC  B»N#MO#M  1»RQ2)  

I F(  A2 1 #NE*  0#  0 ICALL  NINEC N# MO#  P#  PO#  R02#  B#M  1 > 
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IF<M#EQ*NHAX>GOTO  30 


R03-RQ»<  N+ 1) /<  Z 0-<  M* 1 > > 

IF<A21.EQ.0.0)CALL  NBSTXCC#W*  1#M0*M1#R03> 

IFCA2UN  E.0  .0  > CALL  N IN  EC  »♦  UMO#  P*  P0*H03*C#HO 

3 0 1F<HMAX-U0«EQ«M>Q0T0  40  

R04-R0»<  N*  2}/<ZC-<  N+  2 J > _ 

IFCA2UEQ.0.  OKI  ALL  WBSTKC  D#M*  2*M0*M  U R04> 

IFCA21  .NE. 0«0>CALL  NINEC N*2*M0#  P#P0#R04*  D«M1) 

4 0 FN1“C*»B 

FN2-B-A 

FN3-D-C _ 

Q-N 

Q1«B 

Q2-R02  

Q3-N-1 

Q4»A 

Q5-R01 

Q6-N-M 

07«C 

Q8-R03 

I FCFM 1 •GE«0*  0.  AND.  FN2.LE#  0#  O>0OTO  1 00 

IF<N.EQ«KMAX.AND»FN2»LE»0»0)Q0T0  100 

I FCM*  EQ* J#  0# AND* FN 1 • QE« 0*  0> GOTO  100  

O-N-M 

Ql-C 

Q2-R03 


Q4-B 


Q5-R02 

Q6-N+2 

tt7-D 


06»Q7-Q3«0«0 _ 

IF(WWAX-2.0«EQ.M«AND»  FN3.LE*  G«  05  GOTO  1 OC 
Q-l  


IF<N*EQ*2*0*  AMD*  FN2*  GE*  0*  0>GOTO  100. 

XF<FMl*Ea.O>PRINT#« ERROR  »!»■ 

1 F< FM2 * EQ » 0*0)PR1 NT* » ERROR  #2* 

CNl-l/<  l-<  FN3/FN1 ) > 

EN2-1/C  l-<  FN1/FN2>> 
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I  F<  1 ABS<  DM 1 > • QE  • I AB SC  EN  2 > 3N  D-  DM 2 


IF<1ABS<  EN23<QT<2ABS<  EN133ND-EN1 

2F<ND<NE<  03B0T0  98 

IF<FNDGT<Q<  03ND--1 

IF<FNULT<Q<  03ND-2 

98 IF<ND<EQ<DND<<2 

1F<FN1<LT<  0<0<AND<ND<LT<03ND<<  2 

M-M+MD 

1 I F<  N<LE<  0 3N- 1 

IF<M<QT<MMAX3WNMAX 

PR  I NT<N<  B<  R02 

GOTO  5 

1 00  PR2MT200*  Q3<  Q4<  Q5 

PR2NT2  00<  Q<  Q 1 < Q2  

PR2NT200<  Q6<  Q7<  Q3 

2 00  FORMATC  *GW< D20<  1 4<  5X< •NB-*<  D20<  1 4#  /<  • . H01-«<  D20<  1 43 

N-Q 

R0P1-R0 

C0UMT<<C<0 

3 00  ROA-CN* 1 3*R0P1  /<  Z O^N-  1 > 

R0B-N*R0P1/CZ0-N3  

COUNT-COUNT*  1 

CALL  HELPCN+  1<M0<  R0A<AUM1<A<A2 1 3 

CALL  HILP<N<M0<R0B<A3<M  1<B<A2 1 3 

2 F« Aw RQA- BWBOB3<EQ<0<03FH2NT301 

2 FC  < AwR0A-B*R0B3<  EQ<  0<  0<AND<<AUA33<LE<  2E- 143  Q«0<  0 

IF<< A1-A33<GE<-2E- 14<AND<Q<EQ< 0<03ROP1-0<0 
2FCROP1<EQ<0<03Q0T0  310  
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301  FORMATC  * OERROR  #1  •> 

1 FC  COUNT*  EQ*  1 • 0 )QOTO  309 

02— <A1-A33/CCA»R0A-B»B0B3 /ROPJL 3-  DRO  IP 

IFC  02  *LE  • 1 E- 1 4*  AND*  C QUNT  • GE*  5 0.  0 3 ROP  1-0.0 

IFC BOP! *E0* 0*  03G0TQ3 1 0 

3 09  DBQ1P— CA1-A33/CCA»R0A-B»B0B3/R0P13 

IFC COUNT* QT* 2 360T0  307  

IFC DR01P.LT* 0 *0  3C OUNT- 3*0 

IFC C0UNT*EQ*3.0360T0  307 

IFC  COUNT*  EQ.l.  0)  DRO  IP—  DRO  IP 

IFC COUNT. EQ*1 *0)Q7-A1 -A3 

IFC  COUNT *EQ*UO>OOTO  307 

IFCCA1-A3>.GT*0*0*AN D* Q7*LT .0.0 3Q0T0  307 

I FC CA1-A33.LT* 0*0* AND* Q7*0T* 0. 03 GOTO  _3 07 

DRO  IP— DRO  IP 

COUNT- 1*0 

IFC  C ROP 1 ♦DRO 1 P3.LE. 0*  0* AND*  ROP 1 • EQ. 1 E-53GOTO  3 1 0 

3 07  ROP  1 -BOP  1*  DRO  IP 

X FC  ROP  1 *LE*  0*  03 ROP  1 - 1 E-  5 

1 Ft ROP 1 *EQ» IE-5* AND* COUNT* GE*  6.03R0P1-0.0 

IFCRQP1.E0.0.  03G0T0  310 

1 FCDRO 1P*LE*  2E-14* AND* DRO 1 P*  GE*-2E- 1 43 G0T0„3 10 

PRINT* R0P1# DROIP.N.ZO 

QOTO_  300 

3 1 0 R0P2-R0  

COUNT- 1.0 
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COUNT-COUNT* 1 

PRINT*  R OPS* DR02P* N 

» 

GOTO  320 

330 

PRINT340*BOP1 * R0P2 

3 40 

FORMATC "0**  D20* 14* • «R0<  •# D20* 14) 

PRINT900 

9 00 

FORMATC 

1JD0Q_ 

CONTINUE 

STOP 

END 


SUBROUTINE  NIN  EC NTOT* MO*  P»  PO*  RO*NE*M  1 > 

DOUBLE  PRECISION  P<101  >*NTOT,PO*MO» RO*A*C*B* D* E* F* G*H* 2* Y* NB* M 1 

DOUBLE  PREC I SI ON  EL  OG*  DEX  P 

1TEST-NTOT-MO*  1 

IU-ITEST 

I PCI  11  •LT«  1 ) 2 1 1-1  

DO  10  J1-1*I TEST 

J-Jl-l __  

Y-Z-0.0 

IFClTSST*LT«OGOTO  21 

DO  20  Nl-UITEST 

NRESUP— W 1-1 

A-CNBESUP-J>«ELOGCMQ»B 0> 

C-D-0«0 

JFC J.LE. 1 >60T0  35  

JD0_30  K-2*  J 

B-K 
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*•'  • 
K 

*4  ’ , 


I C-C*DLOa<B>  

i IFCNRESUP.LE*  OQOTO  45 

DO  40  K-£«NBESUP 

B-K 

4 0 D-D*DL06<B> 

4 5 E-A*C-D 

CALL  FLOW  E#  J#  IX#  2 Y> 

IF  <E.LT»-174«9963 GOTO  20 

Z -Z  ♦ DEX  PCE  > 

2 0 CONTINUE 

2 1 IF<ITEST»LE.  1>ITEST-1 

ITOT-NTOT^l 

I TEST- 1 TEST* 1 

DO  50  N 1 - I TEST, X TOT 

NBESUP-N1-1 

c-d-o*o 

M-MO 

1 F<  M»LE«  1 >QOTO  65 

DO  60  X-2«M  

B-K 

6 0 c-c*a.oe<B)  

6 5 2 F< J#LE#  1 >80T0  75 

DO  70  K-2jJ 

B-K 

70  I*D*ELO«<B>  

7 5 A-F-0.  0 ___ 

I F<NRESUP«LE*  1 >00T0  35 


A-60 


1 


DO  80  K-2#NRESUP 

B-K 

8  0 A- A* EL OGC  B ) 

8 S M-NTOT-NRESUP 

IFCN#LE#1>6QT0  95 

DO  90  K-2#N 

B-K 

9 0 F-F*EL0e<B> 

9 5 q-<NT0T-M0-J>»EL0GCM0 3 

H-<  NRESUP-J  5* ELOGC  RO)  

E-C+D+G+H-A-F 

CALL  FLOVCE#  J# IX# IY> 

1F<E.LT.-1 74. 996>B0T0  50  

Y-Y4DEXPC  ED 

5 0 CONTINUE 

I F(  «I*EQ*  0)PCN  1 • 0/C  Y*Z  > 

IFCJ.  EQ *0> GO TO  10 

PCJ>1*<CY«-Z) _ 

10  CONTINUE 

CALL  TENCNTOT#MO#P#RO> 

CALL  SIXTENC N TOT# M 1 # F#  NB#  PO  ) 

RETURN  

END 

SUBROUTINE  TENC  NTOT#  MO#  P# RO  > 

DOUBLE  PRECISION  PC  101  >#NT0T#P0#M0#R0#A#B#C# D#E#  F#G#H#Z#  Y 

DOUBLE  PRECISION  CL0Q#DEXP 

INITAL-NTOT-MOM 
I F<  INI TAL«LT*  1 > INI  TAL- 1 


A-  6 1 


ITOT-NTOT 

DO  10  J— IN2  TAL#  1 TOT 


Y-Z-0.0 


DO  20  N1-U1M2TAL 


NRESUP-W1-1 


A-C-H-G-O*  0 


2 FC  J«LE«  1 >GOTO  105 


DO  100  K-2*J 
B-K 


100  0-6*ELOG<B> 


1 05  N-NTOT-d 

IF<N*LE. 1)80T0  115 

DO  110  K-2»N 

B-K 


1 1 0 H-H+ CL06C  B > 

115  I F<NRESUP*LE*  1 XsOTO  35 

DO  30  K-2#NRESUP 

B-K 

3 0 A-A-»d.OQ<  B? 

3 5 M-MO 

IF< ltll» I JGOTO  45 

DO  40  K— 2#N  


B-K 


4 0 
45 


c-oeiog<b> 


D-<NRESUP-J>*ELOG<  RO) 


E-<NRESUP-<NTOT-MO>>»ELOG<MO> 

F-D»E*G*H-A-C 


A-  62 


* 


CALL  FLOW< F»  J# IX# I Y) 

IF  < F*LTf* 1 74.996>60T0  20 

Y-Y^DEXP<  F> 

2 0 CONTI  HUE 

DO  50  NRESUP-IN2 TAL# I TOT 

A— C— D— E— F— 0«  0 

IF< J«LE. 1 JQ0T065 

_ DO  60  K-&»J 

B-K 

6  0 A-A+CLOe<  B > 

6 5 I-WT0T-J  

I F<  I.LE.1)Q0T0  75 

DO  70  K-2»I 

B-K 

7 0 C-OCLO«3<B> 

7 5 I -NTOT-NRESUP 

IFC I.LE. 1 >60 TO  35 

DO  80  K-2«  I 

B-K 

3 0 E-E>  CLOCK  B> 

8 5 IFCNRESUP.LE. 1 >00 TO  95 

DO  90  X-2«NRESUP 
B-K 

9 0 F-F»EL06<B> 

9 5 D-<WBESUP-  J>»CL0Q<  R0> 

B-D+A+C-E-F 

CALL  FLOVC B» J» IX# I Y) 


A-63 


IF<B*LT***174*996)©0T0  50 
Z«Z*DEXP<B) 

5 0 CONTINUE  

P<J5-1.  0/<Y+Z) 

1 0 CONTINUE 

RETURN 

END 

SUBROUTINE  SXXTEM<NT0T#M1#P#MB#P0) 
DOUBLE  PRECISION  NTOT# Ml#  P<  101)#NB#P0#B 

ITOTiNTOT 

NB-O*  

INITAL-NTOT-M 1+ 1 

IFCI N2TAL.LE.03NB-PO 

_ IFCIN1 TAL .LE*  0) INI TAL- 1 

DO  10  NRESUP-INI TAL# I TOT 

B-<  M 1 -<  NTOT-NRSSUP) )»  P< NRESUP) 

1 0 NB-NB*B 

RETURN  

_END 

SUBROUTINE  FLOWCX#J# IX# IY) 

DOUBLE  PRECISION  X 

I F<X>LT»- 1 74*99  6)1 X- 1 X+  1 

IF<X*  QT.  1 7 4 .99  6)1 Y«IY»t __ 

1F<X.LT*-174.996)PRINT5#J  

5 FORMATC  W***UNDERFLOW**»*«  J-*#I3> 


IFCX.QT* 174*996)PRINT10# J 


1 0 FQRMATC  * *****  OVER  FLO  J-a*I3) 

IF<X.QT. 174.9963X-17A.99 

RETURN 

END _ 

SUBROUTINE  NBSTKCNB*NTOT*MO*M  1,  R03 

DOUBLE  PRECI SION  NTOT*Ml*NB* STOCK* D*C*MO* RO# PO 

.DOUBLE..  PRECI  SION  CLOG*  DEXP 

M-NTOT-Ml  + 1 

I TEST*  1 00 

MB— 0*0 

IF<M.LE«0)NB-  DEXP<-MO*RO) 

IF<M«LE»03M- 1 

DO  50  NRESUP-M* I TEST  _ 

C-0.0 

I F<NRESUP*LE#  1 3Q0T0  65 

DO  60  J-2*NRESUP  

D-J 

60  C-C+ELOQCD)  

65 IFCNRESUP*ELO8<MO*RO>-C-MO»RO.LT.-100>eOTO  1 1 1 

5 0 NB«NB+<M1-NT0T+NRESUP)»DEX PC  NRESUP*  DLOGC MO* R_Q3-_C  3* DEXPC  -MO»EO) 

I 1 1 CONTINUE  

RETURN 

END 

SUBROUTINE  HELPC  N*MO*  ROB*  A3*M  1*  B*  A2 1 3 

DOUBLE  PRECISION  N* MO* ROB*  A3*M1*B*A2*F(  1013*  PO*  DERI* A 21  

IFC A21  •EQ*0*03CALL  NBSTKC A2*N-  1*M0*M  1* R0B3  

I F<  A2 1 • EQ« 0* 03CALL  NBSTKC  A3»N*M0*M1*  ROB) 


A-65 


I 


r 

_ _ IFCA2UEa*0.0)DERl-M0»R0B  . . . 

IF<  A21  • ME*  0*  0 ) CALL  MI  NEC  N»  1 # M O#  P*  PO*  ROB*  A£#  M I J> 

I F< A2 1_*  N. E*  0*  0 5 CALL  NINECN#MO#  P*  P0#R0B#A3»M  1 ) 

’ m A21  *NE*0* 0 5 CALL  SUB21CM*P* DERI) 

B--<  DERI/  ROB  >»<  A3-A2> 

RETURN 

.END 

SUBROUTINE  SUB21CN*  P» DERI ) 

DOUBLE  PRECISION  N#P<  1015*DERI  __  

DOUBLE  PRECISION  CLOG#  DEXP 

DERI -CUG 

ITOT-N j 

DO  10  NRESUP-UXTOT 

... 

10  DER1«DER1-*NRESUP*P<NRESUP> 


I 
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impact  of  inventory  investment  level  and  resupply  time  level  upon  backorder 
level.  Specific  solution  procedures  are  developed  for  and  are  applied  to  the 
trade-off  analyses  mentioned  above. 
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